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ABSTRACT 


This  thesis  presents  a  design  criterion  for  calculating 
non-minimum  time  deadbeat  digital  controllers,  for  a  unit 
step  input,  which  minimize  the  noise  transfer  of  a  sampled 
data  system  when  white  noise  is  present.  The  system  consisted 
of  a  linear  plant,  a  digital  controller,  a  zero-order  hold 
and  a  unity  feedback  branch. 

The  non-minimum  time  deadbeat  controllers  are  used  in 
place  of  the  minimum  time  deadbeat  controllers  to  improve 
the  noise  performance  of  the  system.  All  deadbeat 
requirements  are  met  by  the  non-minimum  time  deadbeat 
controllers  except  that  the  error  for  a  unit  step  input  is 
not  reduced  to  zero  in  minimum  time. 

The  design  criterion  is  practical  in  application  as 
it  depends  only  on  the  plant  parameters  and  system  sampling 
rate.  No  explicit  reference  to  the  noise  input  is  contained 
in  the  criterion. 

The  non-minimum  time  controllers  may  have  any  number 
of  extra  terms  as  the  noise  reduction  increases  with  each 
extra  term.  Beyond  a  certain  number  of  extra  terms,  the 
noise  reduction  becomes  uneconomical  because  of  the  extra 
controller  complexity  and  slower  system  response  time.  The 
length  of  this  extension  varies  depending  on  the  design 
requirements  of  the  system.  In  this  thesis,  controllers  with 
up  to  three  extra  terms  were  considered. 
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The  non-minimum  time  controllers  were  tested  using 


both  a  hybrid  and  a  purely  digital  simulation  to 
minimum  noise  transfer  requirement. 
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INTRODUCTION 


When  sampled  data  control  systems  containing  a 
linear  plant  with  a  digital  controller  to  give  deadbeat 
response  to  a  unit  step  input  are  subjected  to  any  noise 
input,  the  errors  which  occur  are  rather  large.  This  is 
caused  by  the  high  gains  required  in  the  digital  controller 
to  satisfy  the  deadbeat  requirements  for  a  unit  step  input. 

These  requirements  are: 

1.  The  system  reaches  zero  steady-state  error  after 
a  minimum  number  of  sampling  periods  with  no 
intersampling  ripple. 

2.  The  system  normally  should  have  no  transient 
overshoots.  In  some  cases  this  overshoot  is 
unavoidable.  Then  the  transient  response  should 
be  as  fast  as  possible,  with  a  finite  settling 
time  measured  at  the  sampling  instants. 

The  controller  which  satisfies  these  requirements 
for  a  given  system  is  unique.  Therefore  if  any  alterations 
are  made  to  the  controller  to  improve  the  noise  performance 
of  the  system,  the  resulting  controller  cannot  satisfy  all 
deadbeat  requirements. 

In  this  thesis,  a  design  criterion  is  presented  which 
yields  controllers  that  minimize  the  noise  transfer  of  the 
system  in  the  presence  of  noise  which  is  assumed  to  be  white. 
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These  controllers  are  designed  to  satisfy  all  deadbeat 
requirements  except  the  minimum  response  time  restriction. 
Therefore  the  predicted  controllers  may  be  called  "non¬ 
minimum  time  deadbeat  digital  controllers." 

When  using  this  method,  the  plant  parameters  and  the 
sampling  rate  of  the  system  are  normally  held  constant. 

Thus  the  non-minimum  time  deadbeat  controller  may  be  directly 
substituted  for  the  corresponding  minimum  time  deadbeat 
controller  for  a  given  system  to  achieve  the  maximum  noise 
reduction . 

This  method  of  extending  the  controller  to  improve 
the  noise  performance  of  the  system  has  the  following 
advantages : 


1.  The  gains  of  the  controller  are  lowered  and  thus 
the  level  of  the  transient  response  is  lowered. 

2.  Because  the  controller  is  no  longer  a  minimum 
time  controller,  extra  degrees  of  freedom  are 
available  in  its  design.  These  may  be  used  to 
minimize  the  transfer  of  noise  by  the  system. 

In  effect,  the  predicted  non-minimum  time  controllers 
are  designed  to  do  as  little  as  possible  to  the  signal 
passing  through  them.  This  is  contrary  to  the  more  standard 
approach.  Usually  the  system  loop  gain  is  made  as  large  as 
possible  to  allow  very  rapid  correction  of  any  disturbances. 
This  is  generally  accomplished  by  the  insertion  of  a  network 
in  the  feedback  path  to  raise  the  loop  gain  without  altering 


•  fcs  i  :o  P.3'1  if laiauBai  ©rto  5o  Iivs  <  : 
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the  forward  path. 

It  would  be  possible  to  insert  a  feedback  compensator 
in  the  deadbeat  system  to  improve  the  noise  performance. 

However  it  is  a  more  complicated  problem  in  terms  of  the 
system  realization  then  is  the  simple  addition  of  terms  to 
the  controller. 
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Authors  such  as  Tou  '  '  # Ackermann  and  Chang  have 

published  methods  for  calculating  controllers  which  minimize 
the  noise  transfer  of  a  sampled  data  system.  However  these 
methods  do  not  readily  lend  themselves  to  the  case  when  deadbeat 
is  an  additional  constraint.  Some  of  the  reasons  for  this 
observation  are  discussed  in  Chapter  3. 

To  the  best  of  the  author's  knowledge,  no  published 
literature  contains  any  specific  reference  to  the  noise 
reduction  problem  for  deadbeat  systems.  Therefore  the  design 
criterion  presented  in  this  thesis  is  derived  independently 
of  published  methods. 

The  proposed  criterion  is  stated  in  Chapter  3.  This 
criterion  may  be  used  to  calculate  controllers  which  have 
any  number  of  extra  terms.  However  at  some  point,  the 
additional  noise  reduction  becomes  uneconomic  when  compared 
to  the  additional  controller  complexity  required. 

The  criterion  presented  is  practical  in  application 
and  it  is  expressed  only  in  terms  of  the  coefficients  of  the 
non-minimum  time  controller.  These  controller  coefficients 
are  functions  of  the  plant  parameters  and  the  sampling  rate 
of  the  system.  The  noise  input  does  not  appear  explicitly 


. 
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in  the  criterion.  To  allow  the  design  criterion  to  take 
this  form,  pure  white  noise  must  be  assumed  as  it  has  no 
cross-correlation  with  a  unit  step  input.  Therefore  the  two 
signals  may  be  assumed  to  be  totally  independent. 

In  Chapters  1  and  2,  the  form  of  the  system 
considered  is  presented  along  with  other  general  information. 
The  details  of  a  hybrid  computer  simulation  used  for  the 
initial  tests  of  the  non-minimum  time  controllers  are 
presented  in  Chapter  3.  In  Chapter  4,  the  details  of  a 
purely  digital  test  simulation  used  to  obtain  the  results 
presented  in  Chapters  5,  6  and  7  are  given. 

The  results  presented  in  Chapters  5,  6  and  7  are  for 
non-minimum  time  controllers  with  one,  two  or  three  extra 
terms  respectively.  In  Chapter  5,  there  is  no  noise  input  to 
the  system  when  the  ramp  and  acceleration  inputs  are  used. 

In  Chapter  7,  results  are  presented  for  the  case 
when  the  sampling  rate  of  the  system  is  not  held  constant. 

The  sampling  rates  used  are  such  that  the  total  response 
times  to  deadbeat  for  both  the  minimum  time  controllers  and 
the  predicted  non-minimum  time  controllers  with  one  extra 
term  are  the  same.  This  approach  overcomes  the  disadvantage 
of  the  slow  response  of  the  non-minimum  time  controller  while 
still  providing  some  noise  reduction. 


. 


CHAPTER  I 


GENERAL  SYSTEM  CONSIDERATIONS 

Throughout  this  thesis,  some  basic  information  on 
the  systems  tested  will  be  required.  For  each  plant  to  be 
considered,  the  following  information  is  presented  in  this 
chapter: 


1.  The  transfer  function  of  the  plant. 

2.  The  -  transform  of  the  plant. 

3.  The  discrete  signal  flow  graph  for  the  system. 

4.  The  discrete  state  equations  of  the  system. 


1-1  SYSTEM  CONFIGURATION 


The  basic  system  considered  is  a  unity  feedback 
system  with  a  noise  contaminated  input  or  feedback  signal 
as  shown  in  FIG.  1.1. 

In  this  system,  the  noise  is  shown  as  injected  in 
the  feedback  path.  Since  this  is  a  unity  feedback  system, 
the  noise  could  actually  be  contaminating  the  input  signal 
or  be  introduced  by  the  feedback  branch. 

The  general  form  of  the  digital  controller  is: 


D  (z) 


h0  +  hl  Z 1  +  h2 


-2  ,  "n 
z  4*  •  *  •  +  hn  z 


(1.1) 


,  ,  -1  ,  -2  ,  -n 

1  +  e±  z  +  e2  z  +  . . .  +  en  z 
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where  n  is  the  total  number  of  sampling  periods  required  to 
reach  deadbeat  for  a  unit  step  input.  The  zero  order  hold 
has  the  transfer  function: 


Gz  (s) 


e 

s 


-st 


(1.2) 


where  T  is  the  sampling  rate  of  the  system.  For  most 
calculations,  the  zero  order  hold  and  the  plant  are  combined 
as  follows: 

G  (s )  =  G7  (s)  G  (s )  (1.3) 

pm  z  m 

where  G^(s)  is  the  particular  plant  under  consideration. 


1-2  TYPES  OF  PLANTS  CONSIDERED 


There  are  several  basic  plant  configurations  which 
represent  the  majority  of  system  classes  possible.  Four 
such  plants  are  considered.  These  are: 

(s )  =  1/s  (s  +  «)  (1.4) 

G2(s)  =  1/s  (s  +  «)  (s  +  3)  (1.5) 

G3(s)  =  (s  +  «)  /s  (s  +  3)  (1.6) 

G4(s)  =  1/s  (s2  +  2?  un  s  +  a)2) 


(1.7) 


'  .  .  '  -o  '. 
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D  (z)  G  pm  (s) 


0 z ( s )  =  ZERO  ORDER  HOLD 

Gm(s)  =  PLANT  UNDER  CONSIDERATION 


FIGURE  1.1  BASIC  SYSTEM  CONFIGURATION  WITH  A  UNIT 

STEP  INPUT  PLUS  NOISE 


(a)  (b) 


FIGURE  1.2 


DEADBEAT  RESPONSE  TO  A  UNIT  STEP  INPUT 
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The  deadbeat  responses  exhibited  by  G^(s)/  an<^ 

G^(s),  whether  for  the  minimum  time  or  extended  controllers, 

are  always  as  shown  in  FIG.  1.2  (a).  However  for  G^(s), 

it  is  possible  for  some  combinations  of  parameters  to  yield 

3 

a  response  as  shown  in  FIG.  1.2  (b) .  This  response  still 

has  no  intersample  ripple  and  reaches  unity  at  nT  as 
specified. 


1-3  TRANSFORMS  OF  PLANTS  CONSIDERED 

As  shown  in  (1.3),  the  plant  Gm(s)  and  the  zero-order 

hold  G  (s)  are  combined  as  G__(s).  For  later  use,  it  is 
z  Pm 

necessary  to  calculate  G  (z) .  This  section  contains  the 
calculation  of  Gp-^(z).  The  transforms  for  the  other  three 
plants  are  then  listed. 

Substituting  (1.2)  and  (1.4)  into  (1.3)  yields: 


-s  T 
e 

(s  +  «) 


(1.8) 


Taking  the  Transform  of  (1.8)  yields: 

1 

s  (st*) 


Pi 


(z>  =  3 


1  -e 


-ST 


s  (s+°0 


=  (1  -  z 


-1 


>  ? 


1-z 


-1 


oc 


T  z 
.  (z-l)2 


—  OC  T 

(1  -  e  ) z 

- -  oc'rp' 

a (z-l)  (z-e  ) 


(1.9) 


s- 
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From  (1.9),  define  the  following: 


A  = 


—  ccT 


e 


(1.10) 


Substituting  (1.10)  into  (1.9)  yields  the  simplified 
expression : 


(z) 


z  (°cT  +  A-l)  +  (1  -  A  -  °cTA) 

a2  (z  -1)  (z  -  A) 


(1.11) 


To  simplify  (1.11),  define  the  following: 

Rl  =  (<*T  +  A  -1)  /  cc2  (1.12) 

R2  =  (1  -  A  -  ocTA)  /  «2  (1.13) 


Substituting  (1.12)  and  (1.13)  into  (1.11),  yields  the 

transform  of  G  , (s)  as: 

Pi 


,  .  _  (Rl)  z  +  R2 
plUJ  "  (z  -1)  (z  -  A) 


(1.14) 


The  transforms  for  the  remaining  three  plants 

are  derived  in  a  similar  fashion. 

Thus  for  Gp2  '  the  transform  is: 


GP2(Z) 


(R7)  z  +  (R8)  z  +  U 
(z  -1)  (z  -A)  (z  -B) 


(1.15) 


Tx) 


:as  (a)  ^0  lo  mioiensizt 
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where 


A 

= 

—  ocT 
e 

(1.16) 

B 

= 

-3T 

e 

(1.17) 

Rl 

= 

1  /  -3 

(1.18) 

R2 

= 

-  (i  /«2e  +  i 

/«B2) 

(1.19) 

R3 

= 

1  /«2  (6  -  «) 

(1.20) 

R4 

= 

i  /e2  (“  -  b) 

(1.21) 

R7 

= 

T  (Rl)  -  R2  (1 

+  A  +  B)  -  R3  (2  +  B)  - 

R4  (2  +  A) 

(1.22) 

R8 

= 

-  (Rl)  T  (A  +  B) 

+  R2  (A  +  B  +  AB)  + 

R3  (1  +  2B)  + 

R4  (  1  +  2A) 

(1.23) 

U 

= 

Rl  (TAB)  -  R2 

(AB)  -  (R3 )  B  -  (R4)  A 

(1.24) 

For 

the  plant,  G^^s) 

,  the  transform  is: 

G 

p3 

(z) 

(B)  z  +  C 

(1.25) 

(z  -1)  (z 

-A) 

where 


A 


(1.26) 
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B  =  (3  (aT  +  A  -1)  +  «  (1  -A))  /oc2  (1.27) 

C  =  (3  (1  ~  A  -  A«T)  (1  -  A))  /«2  (1.28) 


Finally,  the 


transform  of  G  ^  (s) 


is : 


Gp4(z) 


(Al)  z2  +  (A2 )  z  +  A3 _ 

(z  -1)  <z2  +  (AC8)  z  +  AC9) 


(1.29) 


where 


B  =  a)  (1  -  C2)7 
n 


AC1 

= 

1  /  (A2  +  B2) 

AC  2 

= 

-  2  (A)  (AC1 ) 2 

AC  3 

= 

-  AC  2 

AC  4 

= 

-  AC1  -  2  (A)  AC 2 

AC  5 

= 

T  (AC1 ) 

AET 

= 

EXP  (  -  AT) 

BT 

= 

B  (T) 

AC  6 

= 

AC 2  (AET)  COS (BT) 

AC  7 

_ 

((2  (A2)  ACl  -  1)  AC1  (AET)  SIN(BT))/B 

(1.30) 

(1.31) 

(1.32) 

(1.33) 

(1.34) 

(1.35) 

(1.36) 

(1.37) 

(1.38) 

(1.39) 


(1.40) 
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AC 8  =  -  2  (AET)  COS(BT) 


(1.41) 


AC 9  =  (AET) 


Al  =  AC 5  +  AC 2  (AC8  +1)  +  AC 6  +  AC 7 


(1.42) 

(1.43) 


A2  =  AC 5  (AC8 )  +  AC 2  (AC9  -  AC 8)  +  AC 3  - 


2 (AC 6  +  AC7) 


A3  =  AC 9  (AC5  -  AC2 )  +  AC6  +  AC 7 


(1.44) 

(1.45) 


Refer  to  Table  1.1  for  the  listing  of  G  (s) ,  G  (s) 

m  pm 

and  G  (z)  for  all  cases, 
pm 


1-4  THE  FORM  OF  THE  DIGITAL  CONTROLLER 


The  pulse  transfer  function  D(z)  of  the  digital 
controller  is  of  the  form: 


D(z) 


BE 2  (z) 
BE1  (z) 


+  h. 


1  +  e^  z 


+  . 


+  . 


+  h 


n 


+  e 

n 


(1.46) 


where  BEl(z)  and  BE2(z)  represent  the  -  transforms  of 
the  input  and  output  signals  of  the  controller  as  shown  in 
FIG.  1.1.  Cross-multiplying  in  (1.46)  and  taking  the 

4 

inverse  ^  -  transform  yields  the  simplified  form: 

*  n  *  n 
BE 2  (t)  +  E  (e,  )  BE 2  (t  -  kT)  =  E  (h,  ) 

k=l  k  n  K 

BE1* (1  -kT) 


(1.47a) 
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Rearrange  (1.47a)  as  follows: 


n 


BE 2*  (t)  =  Z  (h,  )  BEl*(t  -  kT) 

k=0 


n 


Z  (e,  )  BE 2  (t  -  kT) 
k=l  K 


(1.47b) 


From  (1.47b),  define  the  following: 


n 


y,  (t)  =  Z  (h,  )  BE1  (t  -  kT) 

k=0 


(1.48) 


n 


y9(t)  =  Z  (e,  )  BE 2  (t  -  kT) 

^  k=l  K 


(1.49) 


Rewriting  (1.47b)  using  (1.48)  and  (1.49)  yields  the 
simplified  equation: 


BE 2  (t)  =  Yj^t)  -  y2  (t) 


(1.50) 


Equation  (1.50)  may  be  realized  by  direct  digital 

4 

programming.  This  realization  is  shown  in  FIG.  1.3.  This 
method  of  realization  will  be  used  in  all  cases  considered. 


It  is  used  as  it  lends  itself  easily  to  programming  on  a 
hybrid  or  digital  computer. 


:  t?ajv70l !  o3  srto  9  niiafc  ,  (d.’fr.I)  fno^l 


3)  L:  3  (  1 
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BE2  *(t) 


FIGURE  1.3  DIRECT  DIGITAL  PROGRAMMING  REALIZATION- 

OF  D ( z ) 
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1-5  STATE  EQUATIONS  AND  FLOW  GRAPHS  FOR 
G1 (s)  =  1/s  (s  +  «) 


In  order  to  carry  out  some  of  the  derivations  and 


system  simulations  outlined  later,  it  is  necessary  to 
derive  the  discrete  state  equations  for  the  given  systems. 


For  the  plant  G^s),  the  entire  derivation  is 


presented.  The  signal  flow  graphs  and  discrete  state 
equations  of  the  remaining  systems  will  then  be  presented 
without  derivation. 

Consider  the  system  where  G^(s)  =  1/s  (s  +  a)  is  the 
plant.  The  signal  flow  graph  which  represents  FIG.  1.1  is 
shown  in  FIG.  1.4.4  Applying  Mason's  Rule  to  FIG.  1.4 
yields  the  simplified  form: 


Xx  (s) 


(1.51) 


(1.52) 


Taking  the  inverse  Laplace  Transform  of  (1.51)  and  (1.52) 
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gives : 

x^t)  =  X1(t0)  +  i  (1  -  e-" (t  ■  t0)  )  X2(t0)  + 

l  ( (t  -  tQ)  -  l  (1  -  e-* (t  "  t0)  ) )  BE2  (tg) 

X2(t)  =  (e"“(t  “  V)  X2(tQ)  +  i 

(1  -  e““  (t  "  t0)  )  BE 2  (t*) 

If  t  =  (k  +1)  T  and  =  kT,  then  the  previous  two  equations 
may  be  rewritten  as: 

—  (XT’ 

X1[(k  +1)T]  =  X-JkT]  +  -  (1  -  e  )  X2[kT]  + 

i  (T  -  i  (1  -  e"“T))  BE2 [kT]  (1.53) 

X2[(k  +1) T]  =  (e_“T)  X2(kT]  +  i  (1  -  e_o:T) 

BE2 [kT]  (1.54) 

Using  (1.10),  the  following  constants  may  be  defined  for 
(1.53)  and  (1.54) : 

B  =  (1  -  A)  /  «  (1.55) 


C 


(T  -  B)  / 


CC 


(1.56) 


. 
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Substituting  (1.10),  (1.55)  and  (1.56)  into  (1.53)  and 

(1.54)  yields  the  final  form  of  the  discrete  state  equations 
for  the  system  when  G^(s)  is  the  plant. 

X  [(k  +1)T]  =  Xx[kT]  +  (B)  X2[kT]  + 

(C)  BE2[kT]  (1.57) 

X2[(k  +1)T]  =  (A)  X2[kT]  +  (B)  BE2[kT]  (1.58) 

1-6  STATE  EQUATIONS  AND  FLOW  GRAPHS  FOR 
G2  (s)  =  1/s  (s  +  «)  (S  +  3) 

Applying  the  method  of  Section  1-5  to  FIG.  1.5 
yields  the  discrete  state  equations  of  the  system  containing 
G2 (s)  as  follows: 


X1[(k  +1)T]  =  Xx[kT]  +  (CONS2 )  X2[kT]  + 

(CONS3 )  X3[kT]  +  (CONSl )  BE2 [kT]  (1.59) 

X2  [  (k  +1 )  T ]  =  (A)  X2[kT]  +  (CONS 7)  X^kT]  + 

(CONS  3 )  BE2[kT]  (1.60) 

X3[(k  +1)T]  =  (B)  X3[kT]  +  (CONS9 )  BE2 [kT ]  (1.61) 

Using  (1.16),  (1.17)  and  (1.18),  the  constants  in  the 


■'  !  :>  ( 2 )  r  >  ?  .  9d  t  x .  H 


+  [Trf]  X  (S2M0D)  r  [TJU  lX 


<1  £  ia  (.  3:  a>;  +  {?*}  <  u  cru 
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previous  three  equations  are  defined  as: 

Y  =  a  +  3  (1.62) 

C0NS1  =  R1  (  (T  -  (Rl)y)  +  R1(«2B  -  32  A)/(«  -3)) 


(1.63) 

CONS  2 

=  (1  -  A)  /  « 

(1.64) 

CONS  3 

=  R1  (1  +  (3A  - 

«B)  /  («  -6)  ) 

(1.65) 

CONS  7 

=  (A  -  B)  /  («  - 

6) 

(1.66) 

CONS  9 

=  (1  -  B)  /  3 

(1.67) 

1-7  STATE  EQUATIONS  AND  FLOW  GRAPHS  FOR 
G3Ts)  =  (s  +  3)  /s  (s  +  «) 

From  FIG.  1.6,  the  discrete  state  equations  of  the 
system  incorporating  G3(s)  are: 

X1[(k  +  1 )  T]  =  X1[kT]  +  (CONSl)  X2[kT]  + 

(C0NS2 )  BE  2 [kT]  (1.68) 

X2[(k  +  1)T]  =  (A)  X2[kT]  +  (CONSl) 

BE2[kT]  (1.69) 

The  output  equation  of  the  system  for  this  case  is: 


*  •  ' 

-  T))  Lft 
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C [ (k+1) T] 

=  (8)  x1[(k+l)T]  +  x2[(k+l)T] 

(1.70) 

Using  (1.26)  ,  the 

constants  of  (1.68)  and  (1.69)  are 

defined  as: 

CONS1  = 

(1  -  A)  /  oc 

(1.71) 

CONS 2  = 

(T  -  CONS1 )  / 

(1.72) 

1-8  STATE  EQUATIONS  AND  FLOW  GRAPHS  FOR 

G.  (s)  =  1/s  (s2  +  2CcjOS  +  (ji)2) 

4   n_ n 

The  discrete  state  equations  of  the  system  in 
FIG.  1.7,  which  contains  G^(s)  are: 


X±[  (k  +  l)T]  =  xx[kT] 

+  (CONS1 ) 

X2 [kT] 

+  (CONSC) 

X3[kT] 

+  (CONS2) 

BE2 [kT] 

(1.73) 

X2[(k+1)T]  =  (CONS3 ) 

X2[kT]  + 

(CONSB) 

X3 [kT]  + 

(CONSC ) 

BE 2 [kT] 

(1.74) 

X3[(k+1)T]  =  (CONS  4 ) 

X2[kT]  + 

(CONS 5) 

X3[kT]  + 

(CONSB) 

BE 2 [kT] 

(1.75) 

Using  the  constants  defined  in  (1.30),  (1.31),  (1.32), 


t 

=  S8VCD 
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(1.37)  and  (1.38),  the  constants  in  equations  (1.73), 

(1.74)  and  (1.75)  are  defined  as: 

<(>  =  tan"1  (-  A/B)  *  (1.76) 


CONSA 

= 

(AET) 

COS (BT) 

(1.77) 

CONSB 

= 

(  (AET) 

SIN (BT) )  /B 

(1.78) 

CONSC 

= 

AC1  (1 

-((AET)  COS  (BT  +  <J>)  )/COS  (<J>)  ) 

(1.79) 

CONS1 

— 

CONSB 

+  2  (A)  CONSC 

(1.80) 

CONS  2 

= 

AC1  (T 

+  2  (A)  AC1  (CONSA  -1)  + 

CONSB 

(2 (A2)  AC1  -1) ) 

(1.81) 

CONS  3 

= 

CONSA 

+  (A)  CONSB 

(1.82) 

CONS  4 

= 

2 

-  0) 

n 

CONSB 

(1.83) 

CONS  5 

= 

CONSA 

-  (A)  CONSB 

(1.84) 
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FIGURE  1.4  FLOW  GRAPH  FOR  SYSTEM  CONTAINING  G^s) 


FIGURE  1.5  FLOW  GRAPH  FOR  SYSTEM  CONTAINING  G2(s) 
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FIGURE  1.6  FLOW  GRAPH  FOR  SYSTEM  CONTAINING  G3(s) 


FIGURE  1.7  FLOW  GRAPH  FOR  SYSTEM  CONTAINING  G4(s) 


CHAPTER  II 


HYBRID  COMPUTER  SIMULATION 

In  this  chapter,  the  method  of  simulation,  the 
calculation  procedure  for  the  extended  controllers,  and 
the  testing  procedure  followed  are  detailed. 

When  the  problem  of  noise  reduction  for  the 
deadbeat  controllers  was  first  considered,  there  was  no 
obvious  way  to  predict  the  minimum  noise  transfer  extended 
controllers.  It  was  also  not  known  whether  any  significant 
noise  transfer  reduction  was  obtainable  by  this  extension 
method. 

By  observing  sequential  sets  of  extended  controllers 
in  the  presence  of  noise,  these  two  questions  could  be 
answered.  The  first  method  used  to  perform  these  obser¬ 
vations  was  a  hybrid  computer  simulation.  A  hybrid 
simulation  had  the  advantage  of  using  available  equipment 
and  straight  forward  techniques  to  allow  fairly  rapid  testing 
of  a  variety  of  controllers. 

For  this  simulation,  the  plant  G^(s)  was  considered 
with  one  and  two  extra  terms.  (s)  was  only  considered 

with  one  extra  term. 

By  observing  the  results  for  these  plants,  infor¬ 
mation  which  would  lead  to  a  criterion  for  predicting  the 


minimum  noise  transfer  extended  controllers  was  obtained. 


' 
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2-1  HYBRID  COMPUTER  SYSTEM 

The  computers  used  in  the  experimental  simulation 

were : 


1.  A  Digital  Equipment  Corporation  PDP-8  programmed 
in  FOCAL. 

2.  An  Electronic  Associates  Incorporated  TR-48 
analog  computer. 

3.  A  Redcor  Corporation  Series  610  Linkage  System 
interface . 

The  PDP-8  digital  computer  is  a  4K  core  memory  version  with 
external  magnetic  tape  drives. 

FOCAL,  which  is  an  interpretive  language,  was 
chosen  as  the  programming  language  because  of  the  ease  of 
programming.  Although  FOCAL  is  much  slower  than  machine 
language,  it  was  felt  that  the  ease  with  which  the  new 
values  for  various  controllers  could  be  inserted  in  the 
program  compensated  for  the  slowness.  The  maximum  sampling 
rate  the  program  could  operate  at  was  T  =  0.5  second. 

This  was  not  a  disadvantage  as  the  normal  sampling  rate 
used  was  1  second. 

The  TR-48  was  used  since  it  has  external  mode 
control  on  each  integrator  pair.  This  allows  the  PDP-8  to 
control  the  operating  state  of  any  part  of  the  simulated 
system. 
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2-2  HYBRID  SYSTEM  SIMULATION 


In  order  to  implement  the  problem  on  the  hybrid 
computer,  it  was  first  necessary  to  write  a  digital  computer 
program  to  simulate  the  D(z).  This  digital  program  was  also 
written  to  include: 

1.  Control  of  the  start  and  duration  of  the 
experimental  runs. 

2.  Control  of  the  operating  state  of  the  TR-48. 

3.  Control  of  the  actual  sampling  rate  T. 

4.  The  zero-order  hold  function  which  is  shown 
separately  in  FIG.  1.1  as  Gz(s). 

5.  The  calculation  of  the  root-mean- square  error 
for  the  system  upon  completion  of  the  experi¬ 
mental  run. 

The  derivation  and  listing  of  this  program  is  given  in 
Appendix  A. 

The  measuring  criterion  used  is  the  root- mean-square 
error  of  the  system.  Therefore  it  is  necessary  to  measure 
the  difference  between  the  response  to  a  unit  step  with  and 
without  the  noise  input.  For  a  particular  controller  under 
test,  this  may  be  accomplished  as  shown  in  FIG.  2.1.  In 
FIG.  2.1,  cm(t)  is  the  desired  system  response  to  a  unit 
step  as  shown  in  FIG.  1.2. 

The  test  system  shown  in  FIG.  2.1  is  difficult  to 
implement  as  it  requires  duplication  of  the  system  to  yield 


isiXqmooos  9d  ybib  airid 
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c'  (t) .  It  is  also  somewhat  redundant  as  the  value  of  c  (t) 

m 

is  unity  for  all  but  the  first  n  sampling  periods.  Therefore 
it  is  possible  to  alter  the  hybrid  realization  to  that  shown 
in  FIG.  2.2.  This  yields  the  same  results  as  the  system 
in  FIG.  2.1. 

The  root- mean- square  error  of  the  system  is 

r  i  i 

- 2  7 

[  o'  (t) ]  .  This  assumes  that  the  noise  source  is 

the  same  for  any  controller  measured.  Since  the  noise  source 
may  fluctuate  slightly,  the  root-mean- square  error  is 
redefined  for  this  test.  Using  the  new  definition,  the 
root- mean- square  error  is  actually  the  ratio  of  the  root- 
mean-squares  of  the  output  of  the  system  and  the  noise  input. 
This  allows  direct  comparison  of  the  results  obtained  for 
any  controller. 

Redefining  the  root- mean- square  error  yields: 


ROOT- MEAN- SQUARE  ERROR  = 


[  c'  (t)]2  / 


In  (2.1)  ,  I  c '  <t )  ]  2  and 
[  c'  (t)  3 2  = 

T 


RMSE  = 


[  n (t)  ]  2 


[  n (t)  ]  2 


are  defined  as: 


T 

/0T  [  c'  (t) ]2  dt 


(2.1) 


(2.2) 


>  5o  ot/Iav  srfi  sc  :t neb  nubs:!  Jds  wnawB  oa!&  ax  dH 
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c'(0  = 


FIGURE  2.1  TEST  SYSTEM  TO  MEASURE  SYSTEM  ERROR 

WHEN  NOISE  PRESENT 


C'(r) 


FIGURE  2.2  TEST  SYSTEM  IMPLEMENTED  ON  HYBRID  TO 

MEASURE  SYSTEM  ERROR  WHEN  NOISE  PRESENT 


(  D (Z)  &  ZERO  ORDER  HOLD) 
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FIGURE  2.3  BLOCK  DIAGRAM  OF  HYBRID  COMPUTER  SIMULATION 
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2  1  2 
n  (t)  =  i-  /  1  n2(t)  dt 

T 


(2.3) 


where  TT  =  total  time  of  the  experimental  run.  Substituting 
(2.2)  and  (2.3)  into  (2.1)  yields  in  simplified  form: 


RMSE 


[  o' (t)  ]2 


dt  /  / 


0 


n^  (t)  dt 


(2.4) 


In  the  hybrid  simulation,  the  squaring  and  integration 
are  performed  on  the  analog  computer  using  the  electronic 
multipliers.  At  the  end  of  the  experimental  run,  the  two 
mean-squares  are  A/D  converted  and  the  division  and  square 
root  are  calculated  digitally. 


2-3  STATE  VARIABLE  METHOD  FOR  CALCULATING  MINIMUM 
AND  EXTRA  TERM  CONTROLLERS 

In  Chapter  1,  the  discrete  state  equations  for  the 
four  systems  were  developed.  A  method  for  calculating  the 
extra  term  controllers  is  based  on  these  equations. 

It  assumed  for  these  derivations  that  all  state 
variables  have  a  zero  initial  condition.  That  is: 

X^(0)  =0;  i  =  1,  2,...,m  (2.5) 

The  input  signals  for  these  calculations  are: 


r  (t ) 


u  (t) 


1.0 


(2.6) 


n  (t) 


0 
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(2.7) 


If  the  system  is  to  reach  deadbeat  in  n  sampling  periods, 
then  the  following  values  of  the  state  variables  must  be 
reached. 


X1[  (n 

+  1)T] 

=  1.0 

(2.8) 

X.  [(n 

l 

+  1)  T  ] 

=  0;  i  =2,3,  ...,m 

(2.9) 

The  method  of  calculation  is  derived  for  the  plant 
G1 (s)  =  1/s  (s  +  «)  . 

From  (1.57)  and  (1.58),  the  discrete  state  equations 
for  the  system  containing  G-^s)  are: 

Xx[(k  +1)T]  =  Xx[kT]  +  (B)  X2[kT]  4-  (C) 

BE 2 [kT] 

X2[(k  +1)T]  =  (A)  X2[kT]  +  (B)  BE2  [kT] 

From  FIG.  1.5,  the  quantity  BEl(tJ)  may  be  defined  as: 

BEl  (t  q  )  =  r(tQ)  -  X1(tQ)  =  1  -  X^tp  (2.10) 

As  in  section  1-5,  define  t^  =  kT  to  allow  (2.10)  to  be 
expressed  in  discrete  form.  Thus  (2.10)  becomes: 


BEl [kT] 


1  -  X-^kT] 


(2.11) 
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The  digital  controller  may  be  replaced  in  a 
discrete  state  equation  by  the  equivalent  gain  during  any 
sampling  period. ^  This  equivalent  gain  is  defined  for  the 
kt^1  sampling  period  as: 

K (k )  =  hk  /  ek  (2.12) 

where  hk  and  ek  are  the  controller  coefficients  as  shown  in 
(1.1).  This  equivalent  gain  is  shown  in  FIG.  1.4  to  FIG. 

1.7.  From  these  diagrams  it  may  be  seen  that: 

BE2[kT]  =  K (k)  (BE1 [kT] )  (2.13a) 

Substituting  (2.11)  into  (2.13)  yields  the  controller 
output  in  terms  of  K(k)  as: 

BE2[kT]  =  K  (k)  (1  -  Xx[kT])  (2.13b) 

Substituting  (2.13)  into  (1.57)  and  (1.58)  yields  the 
discrete  state  equations  as: 

X1[(k+1)T]  =  (1  -  (A)  K  (k) )  Xx[kT)  +  (B) 

X2 [kT]  +  (A)  K (k)  (2.14) 

X2[(k  +  1)T]  =  -K(k)  (B)  Xx[kT]  +  (A)  X2(kT]  + 

(B)  K  (k) 

The  minimum  time  controller  for  G^(s)  reaches 


(2.15) 


-  i  V 
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deadbeat  at  k  =  1.  This  controller  is  unique  and  may  be 
calculated  from  (2.14)  and  (2.15)  by  use  of  the  initial 
conditions  in  (2.5),  (2.8)  and  (2.9).  Following  is  the 

calculation  of  this  controller. 

The  state  variables  at  the  end  of  the  sampling 
period  where  k  =  0  are: 


(A)  K ( 0 ) 


(2.16) 


x2[t] 


(B)  K  (  0  ) 


(2.17) 


Similarly  for  k  =  1,  the  state  variables  are: 


X1[2T] 


(1  -  (A)  K(1))X1[T]  +  (B) 


X2[T]  +  (A)  K ( 1 ) 


(2.18) 


X2[2T]  =  -K ( 1 )  (B)  X1[T]  +  (A) 


X2[T]  +  (B)  K ( 1 ) 


(2.19) 


From  (2.8)  and  (2.9),  the  final  values  are: 


X1[2T] 


1.0 


X2  [  2T ] 


0 


Substituting  (2.16),  (2.17)  and  (2.20)  into  (2.18)  and 

(2.19)  yields  the  simplified  state  equations: 


(2.20) 


U)X  (A)  +  [T]  X 
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X1[2T]  =  1.0  =  (A  +  B2)  K ( 0 )  +  (A)  K(l)  + 

(-A2)  K  ( 0 )  K ( 1 )  (2.21) 

X2[2T]  =  0  =  (AB)  K  ( 0 )  +  (B)  K(l)  + 

(-  AB)  K  ( 0 )  K ( 1 )  (2.22) 

Solving  (2.21)  and  (2.22)  yields: 

K  ( 0 )  =  1  /  (A  (1  -A)  +  B2)  (2.23) 

K ( 1)  =  A  /  (A2  -  B2)  (2.24) 

The  expression  (2.11)  defines  the  error  e (t)  shown  in 
FIG.  1.1  when  n(t)  =  0.  From  FIG.  1.1,  this  error  in 

general  form  is: 

e  [T]  =  BE1 [T]  =  1  -  XX[T]  (2.25) 

It  can  be  shown  that  the  denominator  terms  of  a  deadbeat 
digital  controller  are  the  same  as  the  system  error. 

Therefore  e^  of  the  controller  may  be  written  as: 

ex  =  e [T]  =  1  -  X1 [T]  (2.26) 

Substituting  (2.16)  and  (2.23)  into  (2.26)  and  solving 
yields : 

e±  =  A  /  (A  (1  -  A)  +  B2)  (2.27) 


■ 


3 


35 

The  coefficients  of  the  minimum  time  deadbeat  digital 
controller  for  G^s)  =  1/s  (s  +  °0  may  be  calculated  from 
(2.12),  (2.23),  (2.24)  and  (2.27)  to  yield  in  simplified 

form: 

h  0  =  K ( 0 )  =  1  /  (A  (1  -A)  +  B2)  (2.28) 

h±  =  e±  K  ( 1 )  =  A2  /  (A  (1  -  A)  +  B2) 

(A2  -  B2)  (2.29) 

ex  =  A  /  (A  (1  -  A)  +  B2)  (2.30) 

where  A  and  B  are  defined  in  (1.16)  and  (1.17). 

Since  the  minimum  time  deadbeat  controller  is  unique, 
no  noise  reduction  can  be  performed  on  it  without  altering 
its  form.  The  method  considered  involves  adding  extra  terms 
to  the  controller  without  losing  the  other  deadbeat 
characteristics.  In  this  way,  additional  degrees  of 
freedom  in  the  design  of  the  controller  are  available.  These 
may  possibly  be  used  to  improve  the  noise  performance  of  the 
system. 

Because  additional  degrees  of  freedom  are  available 
in  the  controller  design,  an  infinite  number  of  possible 
controllers  satisfy  the  extended  deadbeat  requirements. 

Some  method  of  predicting  which  controllers,  if  any,  have 
the  best  noise  performance  must  be  derived.  This  is  where 
the  hybrid  simulation  is  useful.  The  performance  of  various 


. 
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extended  controllers  may  be  appraised  in  the  presence  of 
noise . 
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The  deadbeat  controller  which  is  extended  by  one 
term  has  the  following  form  for  the  system  containing 
G1(s)  : 


h  +  hn  z  ^  +  h9  z  2 

D  (z )  =  —2 - ^ -  (2.31) 

1  +  e^  z  +  z 


In  (2.31),  anc^  h^  are  the  extra  controller  terms.  The 
following  derivation  will  show  that  it  is  possible  to 
express  h  q,  h^,  h 2  and  e^  as  functions  of  e^,  «  and  T. 

By  extending  the  controller  by  one  extra  term,  one  extra 
degree  of  freedom  is  available  as  e^  becomes  the  only 
unknown . 

For  a  deadbeat  controller,  e^  is  the  system  error 
at  the  end  of  the  first  sampling  period.  Because  of  the 
deadbeat  constraint,  e^  must  always  lie  between  0  and  1.0. 
All  of  the  controllers  represented  by  (2.31)  may  thus  be 
solved  by  inserting  all  values  of  e^  into  (2.28)  to  (2.30). 
This  is  of  course  not  practical. 

Instead  a  small  number  of  controllers  with  various 
consecutive  values  of  e^  will  be  tested  on  the  hybrid. 
Therefore  extended  controllers  of  the  form  (2.31)  will  be 
solved  for  by  "iterating"  e^  over  the  range  0  ^  ei  ^  1*0 
for  given  values  of  a  and  T. 

As  for  the  minimum  time  deadbeat  controller,  the 
discrete  state  equations  are  written  for  the  various 
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sampling  periods.  For  k  =  0,  the  state  variables  are: 

X1[T]  =  (A)  K ( 0 )  (2.32) 

X2[T]  =  (B)  K  ( 0 )  (2.33) 

For  k  =  1,  the  state  variables  are: 

X]_[2T]  =  (1  -  (A)  K(1))X1[T]  +  (B) 

X2[T]  +  (A)  K  (1 )  (2.34) 

X2[2T]  =  -K  (1 )  (B)  XX[T]  +  (A) 

X2[T]  +  (B)  K ( 1 )  (2.35) 

For  k  =  2 ,  the  state  variables  are: 

X1[3T]  =  (1  -  (A)  K ( 2 ) )  X1[2T]  +  (B) 

X2[2T]  +  (A)  K ( 2 )  (2.36) 

X2[3T]  =  -K (2)  (B)  X1[2T]  +  (A) 

X2[2T]  +  (B)  K ( 2 )  (2.37) 

The  final  values  of  the  state  variables  are: 


X^T]  =  1.0 

X2[3T]  = 


0 


(2.38) 


Substituting  (2.32)  to  (2.35)  and  (2.38)  into  (2.36)  and 
(2.37)  yields  in  simplified  form: 
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X1[3T]  =  1.0  =  (A  +  B2  +  AB2)  K(0)  + 

(A  +  B2)  K ( 1 )  +  (A)  K (2 )  +  (-A2  -  AB2) 

K  ( 0 )  K  (1 )  +  (-A2  -  AB2)  K  ( 0 )  K(2)  + 

(-A2)  K  ( 1 )  K  (2  )  +  (A3)  K  ( 0 )  K  ( 1 )  K(2) 

(2.39) 

X2[3T]  =  0  =  (A2B)  K(0)  +  (AB)  K  (1)  +  (B) 

K  (2 )  +  (-A2B)  K(0)  K  (1 )  +  (-AB  -B3) 

K ( 0 )  K ( 2 )  +  (-AB)  K  ( 1 )  K  (2 )  + 

(A2  B)  K  ( 0 )  K  ( 1 )  K  (2  )  (2.40) 


Eliminate  K(2)  and  all  its  products  from  the  above 
expressions  in  the  following  manner: 3 


1. 

MULTIPLY  (2.39)  by  B 

2. 

MULTIPLY  (2.40)  by  -A 

3. 

Adding  the  results  of  (1) 

and  (2)  yields 

the 

simplified 

expression: 

B 

=  K ( 0 )  (A  +B2  +  AB2  -A3) 

(B)  +  K  ( 1 )  (A  + 

B2 

(B)  +  K  ( 0 )  K  ( 1 )  (A2  -  AB2  +  A3)  (B) 


(2.41) 


. 


39 


Solving  for  K(l)  in  terms  of  K(0)  from  (2.41)  yields: 

K ( 1 )  =  (1  -  K ( 0 )  (A  +  B2  +  AB2  -  A3))  / 

( (A  +  B2  -  A2)  +  K ( 0 )  (A3-AB2  -  A2))  (2.42) 

K(2)  may  be  expressed  from  (2.40)  in  terms  of  K(0)  and  K(l) 
in  the  following  form: 

K  ( 2 )  =  (K  ( 0 )  (-  A2B)  +  K  (1)  (-  AB)  +  K(0)  K(l) 

(A2B))  /  (B  +  K ( 0 )  (-  AB  -  B3)  + 

K ( 1 )  (-  AB)  +  K ( 0 )  K (1)  (A2B) )  (2.43) 

As  stated  in  (2.26): 

ek  =  1  -  X1[kT]  (2.44) 

Substituting  for  k  in  (2.44)  yields  the  following 
expressions : 

e±  =  1  -  X1[T]  (2.45) 

e2  =  1  ~  XX[2T]  (2.46) 

Substituting  (2.32)  into  (2.45)  yields  K(0)  as  the  following 
function  of  e^: 


K  ( 0 ) 


(1  -  ex)  /A 


(2.47) 


:imoi  pniwolloi  ni 
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If  is  iterated  over  the  range  0  to  1,  all  possible 
values  of  K(0)  are  obtained.  These  in  turn  allow  the 
calculation  of  all  possible  values  of  K(l)  and  K(2)  from 
(2.42)  and  (2.43).  Substituting  (2.34)  into  (2.46) 
yields  e 2  as  a  function  of  K(0)  to  K(2)  which  are  now 
functions  of  e^.  The  coefficient  e ^  is  obtained  as: 

e2  =  1  +  K  ( 0)  (-A  -B2)  +  K  (1)  (-A)  + 

K  ( 0)  K  (1)  (A2)  !  (2.48) 

Table  2.1(a)  lists  the  coefficients  of  the  minimum 

time  deadbeat  for  G, (s) .  The  extended  controller  for  G, (s) 

-L  r  '  ■  1 

is  given  in  Table  2.1(b).  These  coefficients  in  Table  2.1(b) 
can  all  be  expressed  as  functions  of  e^. 


k 


hk  =  K  (k)  (ek) 


0 


_ 1 _ 

(A  (1  -  A)  +  B2) 


1 


1 


(A  (1-A)+B2)(A2-B2) 


61  = 


(A(l-A)+B  ) 


TABLE  2.1(a)  COEFFICIENTS  OF  MINIMUM  DEADBEAT 
CONTROLLER  FOR  G-^s)  =  1/s  (s  +  «) 


•c 
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The  resultant  equations  listed  in  Table  2.1(b)  were 
programmed  in  Fortran  and  the  appropriate  iteration  in  e^ 
carried  out.  This  calculation  yielded  a  set  of  controllers, 
with  one  extra  term,  for  simulation  on  the  hybrid  system. 


2-4  z-DOMAIN  METHOD  FOR  CALCULATING  MINIMUM 

AND  EXTRA  TERM  CONTROLLERS 

The  state  variable  method,  outlined  in  Section  2-3, 
is  too  cumbersome  for  higher  order  systems  and  thus  another 
method  was  used.  This  method  is  based  on  the  z-domain 
system  transfer  function.  For  a  unity  feedback  system,  the 
system  transfer  function  is  of  the  form: 


C  (z) 
R  (z) 


D ( z )  G  (z) 

pm  _ 

1  +  D(z)  G  (z) 

pm 


(2.49) 


It  can  be  shown  for  a  system  with  deadbeat  response  to  a 
unit  step,  that  the  conditions  on  the  state  variables  given 
in  (2.8)  and  (2.9)  are  met  if  the  following  equation  is 
satisfied  : 


1  +  D(z)  G  ( z )  =  zn  (2.50) 

pm 

where  n  is  the  number  of  sampling  periods  required  to  reach 

5 

deadbeat . 

To  illustrate  the  use  of  this  method,  consider  the 
example  of  G2 (s)  =  1/s  (s  +  «)  (s  +  8)  .  From  (1.15),  the 


.  3UO  b  J/IT-BO 


(1.15) 


transform  of  G^fs)  is: 

r  -  (R7)  z2  +  (R8)  z  +  U 

p2  (z  -1)  (z  -A)  (z  -B) 

Assume  the  minimum  deadbeat  controller  is  the  following 
form: 


hn  (z  -  A)  (z  -  B) 

D  ( z )  =  -  = 

z  +  e^  z  +  e2 

h0  -  ho(A  +B)  Z_1  +  ho(AB)  z’2 

7“  -1  ,  -2  (2.51) 

1  +  e^  z  +  z 

Therefore  the  numerator  coefficients  of  the  controller  are 
defined  from  (2.51)  as: 

h0  =  h0 

h±  =  -hQ  (A  +  B) 

h2  =  hQ (AB)  (2.52) 

The  digital  controller  is  assumed  to  be  the  form  (2.51), 
because  for  deadbeat  to  occur  the  D(z)  must  cancel  all  the 

c 

poles  of  G  (z)  except  the  pole  at  z  =  -1. 

pm 

For  G2(s),  the  minimum  deadbeat  occurs  at  n  =  3 . 
Substituting  (1.15)  and  (2.51)  into  (2.50)  yields  the 


«  . 
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simplified  form: 

z3 [1  +  hQ (R2  +  R3  +  R4) ]  +  z2 [hQ (Rl)T  -  hQ (R2) 

(1  +  A  +  B)  -  hQ  (R3)  (2  +  B)  -  hQ  (R4) 

(2  +  A)  +  (e±  -1)]  +  z[~  h0(Rl)T  (A  +  B)  + 

hQ  (R2 )  (A  +  B  +  AB)  +  hQ  (R3)  (1  +  2B)  + 

hQ  (2A  +  1)  R4  +  (e2  -  e^  ]  +  [hQ  (Rl)  T  (AB)  - 

hQ (R2 )  AB  -  hQ  (R3)B  -  hQ (R4)  A  -  e2 ]  = 

z6  +  [e^  -1]  z5  +  [e2  -  e^]  z4  - 
[e2]  Z3  (2.53) 

Equating  the  coefficients  of  zn  in  (2.53)  results  in  the 
following  relationships: 


For  n  =  3 : 

1  +  hQ  (R2  +  R3  +  R4)  =  0  (2.54) 

For  n  =  2: 

h0(Rl)T  -  h0(R2)  (1  +  A  +  B)  -  h0(R3)  (2  +  B)  - 

h0(R4)  (2  +  A)  +  (ei  -1)  =  0  (2.55) 

For  n  =  1 : 

-h0(Rl)T  (A  +  B)  +  hQ(R2)  (A  +  B  +  AB)  +  hQ (R3) 

(1  +  2B)  +  hQ  ( 2 A  +  1)  R4  + 


0 


(2.56) 


• 
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For  n  =  0 : 

h  (R1)T  (AB)  -  hQ(R2)AB  -  hQ (R3 ) B  -  hQ (R4)A 

-  e2  =  0  (2.57) 

To  simplify  the  manipulation  of  these  equations,  define  the 
following  terms: 


R5  =  (R1)T  (AB)  -  (R2)AB  -  (R3)B  -  (R4)A  (2.58a) 

R6  =  R4  (2  +  A)  +  R3  (2  +  B)  +  R2  (1  +  A  +  B) 

-  (R1)T  (2.58b) 

Solving  (2.54)  to  (2.57)  and  using  (2.58)  to  simplify  the 
result  yields  the  controller  coefficients  as: 

hQ  =  1/(R5  -  R6  -  (Rl )  T  (A  +  B)  +  R2  (A  +  B  +AB) 

+  R3  (2B  +  1)  +  R4  (1  +  2A ) )  (2.59) 


e1  =  1  +  hQ  (R6 ) 


(2.60) 


e2  =  hQ  (R5 )  (2.61) 

These  results  are  tabulated  in  Table  2.2(a). 

The  controller  with  one  extra  term  is  derived  in  a 
similar  manner.  As  in  Section  1-3,  the  value  of  e^  is 


iterated  between  0  and  1  to  solve  the  extra  unknown. 
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Equation  (2.50)  when  the  extended  controller  is  sought 
becomes : 


4 

1  +  D (z)  Gp2  (z)  =  z  (2.62) 

Assume  the  following  controller  form: 

hQ  (z  -  A)  (z  -B)  (z  -  C) 

D  ( z )  =  — ^ - 2 -  = 

z  +  e^  z  +  e2  z  +  e^ 

h0  +  hiz_1  +  h?z~2  +  h3z~3 

— - - “ - ^2 - - - —  (2.63) 

1  +  e^  z  +  e2  z  +  e^  z 

where  C  is  the  additional  unknown  zero  of  D(z). 

Defining  the  numerator  coefficients  of  the  controller  from 
(2.63)  yields: 


h±  =  -  hQ  (A  +  B  +  C) 
h2  =  hQ  (AB  +  AC  +  BC) 


h 


3 


-h  Q  ABC 


(2.64) 
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Substituting  (1.15)  and  (2.63)  into  (2.62)  yields  the 
simplified  result: 

z4  +  z3[e1  +  hQ(R7)  -1]  +  z2[e2  -  e1  +  hQ(R8)  - 

hQ  C (R7 ) ]  +  z[e3  -  e2  +  hQ  (U)  - 

hQ  C(R8) ]  +  [-e3  -  h0CU]  = 

z8  +  z7[e1  -  1]  +  z6[e2  -  e1l  + 

z5[e3  -  e2 ]  +  z4[-e3]  (2.65) 

3 

Equating  the  coefficients  of  z  yields: 

el  +  h0  -1  =  0 

From  this  equation  h^  is  calculated  to  be: 

hQ  =  (1  -  e±)  /  R7  (2.66) 

Since  e^  is  iterated  from  0  to  1 ,  it  is  now  only  necessary 
to  solve  for  the  other  coefficients  in  terms  of  hg.  Equating 
the  coefficients  of  zn  in  (2.65)  and  solving  yields: 


C  =  1  -  1/ (hQ (R7  +  R8  +  U) ) 


(2.67) 


e2  =  hQ  (U  -  C(R8  +  U)  ) 


(2.68) 


e3  =  -hQ  CU 


(2.69) 


,.9-j  ■  I(8H)D  vrl 


0 


1/ (R5-R6- (R1)T (A+B)  + 


1.0 


1 

2 


h 


0 


R2  (A+B+AB)  +  R3(2B+1)  + 
R4 (2A  +1) ) 

h±  =  -hQ (A  +  B) 
h2  =  h0AB 


e1  =  1  +  hQ (R6) 
e2  =  h0(R5) 


TABLE  2.2(a)  COEFFICIENTS  OF  MINIMUM  DEADBEAT  CONTROLLER 

FOR  G2(s)  =  1/s  ( s  +  «)  (S  +  3) 


k 

0 

1 

2 

3 


hQ  =  (1  -  ei)/R7 

hl  =  _h0  <A  +  B  +  O 

h  2  =  hQ (AB  +  AC  +  BC) 

h3  =  -hQ (ABC) 


fk 

1.0 


0  ^  e1  ^  1  (ITERATED) 


e2  =  hQ(U  -  C 


(R8  +  U) ) 


e3  =  _h0  CU 


TABLE  2.2(b)  COEFFICIENTS  OF  4  TERM  DEADBEAT  CONTROLLER 

FOR  G2(s)  =  1/s ( s  +  «) (s  +  3)  (ONE  EXTRA  TERM) 


i  ;JL10  'MOO  TA  a.  I  •v/MM  M  10  3-^2  C  '.'r  iOD 


49 


These  results  are  summarized  in  Table  2.2(b).  All  these 
coefficients  are  calculable  using  e^  as  the  only  extra 
unknown . 

2-5  NOISE  SOURCE  USED  IN  HYBRID  SIMULATION 

The  random  noise  generator  used  to  obtain  the 

experimental  data  is  a  SERVOMEX  TYPE  R.G.  77.  The  noise 

spectrum  of  this  generator  is  shown  in  FIG.  2.4.  The  low 

frequency  output  of  this  generator  is  flat  to  d.c.  according 

to  the  manufacturer's  specifications. 

It  is  known  from  sampling  theory  that  any  continuous 

input  signal  can  be  represented  by  the  bandwidth  -f  /2  to 

s 

f  /2  after  sampling,  where  f  is  the  sampling  frequency, 
s  s 

Also  the  input  signal  must  be  strictly  bandlimited  between 

-f  /2  and  f  /2  for  the  sampling  process  to  uniquely  reconstr- 
s  s 

uct  the  input. ^  If  the  input  signal  is  pure  white  noise 
and  is  strictly  bandlimited,  the  output  of  the  sampler  is 
bandlimited  white  noise.  The  auto-spectral  density  is 
shown  in  FIG.  2.5  for  bandlimited  white  noise.7 

No  suitable  low  pass  filter  was  available  for  this 
experiment.  There  was  however  a  Krohn-Hite  330M  tunable 
bandpass  filter  available.  This  filter  has  a  slope  of  -80 
db/decade  and  a  tunable  range  from  0.2  Hz  to  20  KHz.  In 
order  to  generate  an  approximation  to  low-passed  white 
noise,  the  aliasing  properties  of  the  sampling  process  are 


used . 
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HZ 


FIGURE  2.4  NOISE  SPECTRUM  OF  SERVOMEX  R.G.  77 

RANDOM  NOISE  GENERATOR 


i 

,  AUTO-SPECTRAL  DENSITY 

_ w  f 

-  fs/2  0 

-  fs/2 

FIGURE  2.5  AUTO-SPECTRAL  DENSITY  OF  BAND LIMITED 

WHITE  NOISE  AFTER  SAMPLING 
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If  the  sampling  rate  of  the  system  is  1.0  second, 

the  value  of  1/2  is  0.5  Hz  as  f  =  1.0  Hz.  Therefore  to 

s  s 

keep  the  aliased  frequencies  from  adding  to  the  low  frequencies 
already  present,  the  lower  filter  breakpoint  is  set  at  1.0 
Hz.  The  upper  breakpoint  is  arbitrarily  set  at  10Hz.  Using 
this  passband,  the  input  signals  at  0.5  Hz  are  at  least  24  db 
down  compared  to  the  middle  of  the  band. 

The  Nyquist  frequency  is  defined  as 

fc  =  fs  /2  (2.70) 

For  any  frequency  f  where  0  ^  f  ^  f  ,  the  higher  frequencies 

o 

7 

that  aliase  with  f  are: 

(2f  ±  f ) ,  (4f  ±  f),  ...  ,  (2Nf  ±  f ) ,  ...  (2.71) 

| 

Each  frequency  f  has  an  equal  number  of  higher  frequencies 

aliased  with  itself.  Since  the  bandwidth  of  0  to  0.5  Hz  is 

narrow,  the  contributions  of  the  aliased  frequencies  should 

be  approximately  constant  over  this  range.  Therefore  since 

the  amount  of  input  signal  below  f  is  negligible,  the 

c 

resultant  output  of  the  sampler  should  be  a  usable 
approximation  to  bandlimited  white  noise  between  0  and  0.5  Hz. 


2-6  TEST  PROCEDURE  ON  THE  HYBRID  SIMULATION 

The  noise  generated  by  the  method  outlined  in  Section 
2-5  is  applied  to  the  system  shown  in  FIG.  2.3.  Experimental 
runs  of  100  seconds  duration  were  taken.  A  run  of  100 


io  wpei  wo  I  o*  pnib^B  moi5  e-  i  ys  bui.  it  £  fci*  id}  qs9?( 
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seconds  allows  the  measurement  of  at  least  one  complete  cycle 
of  all  frequencies  above  0.01  Hz. 

Every  additional  decade  below  this  requires  an 
extra  factor  of  10  in  the  length  of  the  run.  The  additional 
accuracy  obtained  could  be  offset  by  difficulties  with  the 
total  drift  and  overall  accuracy  of  the  analog  components. 

The  hybrid  test  procedure  is  as  follows: 

1.  The  minimum  time  deadbeat  controller  is  tested 
with  a  unit  step  to  verify  the  deadbeat  response. 

2.  This  controller  is  then  tested  in  the  presence 
of  noise. 

3.  The  extended  deadbeat  controllers  are  tested 
with  a  unit  step  to  verify  the  deadbeat  response. 

4.  These  extended  deadbeat  controllers  are  tested 
in  the  presence  of  noise. 

Whenever  a  controller  was  tested,  7  to  10  runs  of 
100  seconds  in  length  were  taken.  The  mean  of  these  runs 
was  calculated  along  with  the  95%  confidence  limits  on  the 
mean.  These  results  are  presented  graphically  in  Chapter 
3.  The  interpretation  of  these  results  and  a  proposed 
prediction  criterion  for  the  extended  controllers  are  also 
given  in  Chapter  3.  The  hybrid  system  errors  are  common  to 
both  the  minimum  time  and  extended  controllers  and  therefore 
the  errors  tend  to  cancel  each  other  when  the  ratio  of  the 


noise  transfers  is  taken. 


« 
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CHAPTER  III 


HYBRID  SIMULATION  RESULTS  AND  THE  DERIVATION 

OF  A  DESIGN  CRITERION 

In  this  chapter,  the  results  and  conclusions 
obtained  from  the  hybrid  computer  simulation  described  in 
Chapter  2  are  presented.  The  results  obtained  from  the 
hybrid  simulation  are  presented  graphically  and  then  some 
preliminary  conclusions  are  stated.  Some  details  of  two 
criteria  considered  are  then  presented. 

The  final  section  contains  the  reasoning  behind  and 
the  derivation  of  the  proposed  design  criterion.  The 
controllers  derived  by  use  of  this  criterion  are  then  compared 
with  the  controllers  which  minimize  the  noise  transfer  of  the 
system  when  tested  on  the  hybrid. 


3-1  SYSTEMS  TESTED  BY  HYBRID  SIMULATION 


Using  the  methods  outlined  in  Sections  2-3  and  2-4, 
controllers  with  one  extra  term  were  calculated  by  iterating 
e^.  The  plants  considered  were: 


1 .  ( s ) 

2.  G1(s) 

3.  G1(s) 

4 .  G2  ( s ) 


1/s (s  +1)  with 
1/s (s  +3)  with 
1/s (s  +1 )  with 
1/s (s  +1)  (s+  2) 


T  =  1  second 

T  =  1  second 

T  =0.5  second 
with  T  = 


1  second. 


' 
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The  noise  transfers  measured  for  the  systems 
containing  each  controller  are  plotted  versus  e^  which  is 
the  iterated  variable.  All  noise  transfers  of  the  systems 
using  the  extended  controllers  are  normalized  with  respect 
to  the  noise  transfer  of  the  system  using  the  appropriate 
minimum  time  controller. 

Controllers  extended  by  two  terms  were  also  considered 
for  G^(s)  =  1/s  (s+1)  with  T  =  1  second.  In  this  case, 

e^  and  e ^  were  the  unknown  variables.  To  simplify  the 
problem  of  calculating  these  controllers,  e ^  was  fixed  at 
0.50  and  then  e^  was  iterated.  The  normalized  noise  transfers 
for  this  system  are  then  plotted  versus  e^  with  =  0.50. 

3-2  GRAPHICAL  RESULTS  AND  SUMMARY 

As  stated  in  Chapter  2,  the  mean  noise  transfer  for 
a  given  system  was  the  average  of  7  to  10  test  runs.  The 
95%  confidence  limits  were  also  calculated  for  this  mean  value. 
For  each  minimum  time  controller,  the  mean  and  confidence 
limits  of  the  noise  transfer  for  the  system  are  listed  in 
Table  3.1. 

FIG.  3.1  to  FIG.  3.5  present  the  normalized  noise 
transfers  for  the  systems  containing  the  extended  controllers. 
The  normalized  noise  transfer  was  calculated  by  dividing  the 
mean  noise  transfer  obtained  using  the  extended  controllers 
by  the  mean  obtained  using  the  appropriate  minimum  time 
controller.  The  95%  confidence  limits  were  then  calculated 


around  the  normalized  mean. 


•- 
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This  method  of  presentation  is  used  to  show  more 
clearly  the  actual  improvement  due  to  the  extended  controllers 


Plant 


T 

Seconds 


Mean  Noise 
Transfer 
Measured  On 
Hybrid  (%) 


95% 

Confidence  Limits 
_ On  Mean _ 

Lower  Upper 


s (s  +1) 


1.0 


88.686 


83.359 


94.013 


s  (s  +3) 


1.0 


70.982 


64.853 


77.112 


s  (s  +1) 


0.5 


146.18 


139.47 


152.89 


s  (s  +l)  (s  +2 )  1 . 0 


95.876 


84.398 


107.36 


TABLE  3.1  NOISE  TRANSFERS  AND  CONFIDENCE  LIMITS  FOR 

SYSTEMS  CONTAINING  THE  MINIMUM  TIME  CONTROLLERS 
MEASURED  ON  THE  HYBRID 


A  summary  of  the  results  shown  in  FIG.  3.1  to  FIG. 
3.5  is  given  in  Table  3.2.  The  results  are  given  as  the 
range  of  e^  which  yields  the  minimum  normalized  noise 
transfer  and  the  range  of  the  normalized  noise  transfers 
corresponding  to  the  e^  values. 


<  .0 
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FIGURE  3.1  HYBRID  RESULTS  FOR  3  TERM  CONTROLLERS  OF  G1  ( s)  =  1/s  (s 

WITH  T=1  SEC. 


NORMALIZED  NOISE  TRANSFER 
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©1 


FIGURE  3.2  HYBRID  RESULTS  FOR  3  TERM  CONTROLLERS  OF 

G1(s)  =  1/s  ( s  +  3)  WITH  T=1  SEC. 

HYBRID  RESULTS  FOR  3  TERM  CONTROLLERS  OF 
G1(s)  =  1/s (s  +  1)  WITH  T=0 . 5  SEC. 


FIGURE  3.3 


NORMALIZED  NOISE  TRANSFER 
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FIGURE  3.4  HYBRID  RESULTS  FOR  4  TERM  CONTROLLERS  OF 

G2(s)  =  1/s ( s  +  1) (s  +  2)  WITH  T=1  SEC. 

FIGURE  3.5  HYBRID  RESULTS  FOR  4  TERM  CONTROLLERS  OF 

G  (s)  =  1/s  ( s  +  1)  WITH  T=1  SEC.  AND  e  =0.5 
1  ^ 
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From  the  summary  in  Table  3.2,  a  number  of  preliminary 
conclusions  may  be  drawn: 

1.  It  is  possible  to  minimize  the  noise  transfer 
of  a  system  containing  an  extended  controller 
with  a  particular  value  of  e^. 

2.  The  improvement  is  significant  considering  the 
added  controller  complexity  when  only  one  term 
is  added.  Comparing  the  results  in  FIG.  3.1 

and  FIG.  3.5  does  not  indicate  if  the  improvement 
obtained  by  two  extra  terms  is  significant. 

3.  Generally  the  range  of  values  of  e^  over  which  a 
minimum  may  occur  is  relatively  wide. 

Considering  these  results,  it  seems  worthwhile  to 
study  this  method  of  noise  improvement  in  more  detail. 

3-3  CONSIDERATIONS  INVOLVED  IN  THE  DEVELOPMENT 

OF  A  DESIGN  CRITERION  FOR  AN  EXTENDED  CONTROLLER 

The  results  given  in  Section  3-2  were  obtained 
experimentally.  This  method  of  finding  the  minimum  noise 
transfer  extended  controllers  is  time  consuming.  A  quicker 
method  would  be  a  design  criterion  to  predict  the  minimum 
noise  transfer  extended  controllers. 

As  stated  in  the  introduction,  the  application  of 
the  criterion  to  a  particular  plant  should  be  as  practical 
as  possible.  This  indicates  that  the  criterion  should  be 
relatively  simple  in  form.  The  simplest  form  of  criterion 
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would  be  in  terms  of  the  plant  parameters  with  no  direct 
reference  to  the  noise  input. 

If  the  input  noise  is  white,  this  assumption  is 
possible.  Stationary  white  noise  has  a  constant  auto-spectral 
density  and  may  be  assumed  to  have  no  cross-correlation  with 
a  unit  step.  Providing  that  the  noise  is  at  least  a  reason¬ 
able  approximation  to  white  noise,  it  should  be  possible  to 
neglect  it  in  the  derivation. 

Several  criteria  were  proposed  and  tested.  The 
predicted  results  of  each  case  were  considered  to  verify  the 
validity  of  the  proposal. 

3-4  CRITERIA  CONSIDERED 

The  first  method  considered  did  not  exclude  the 
noise  input  as  suggested  in  Section  3-3.  This  approach  was 
based  on  the  work  by  Tou^"  and  corrected  by  Jury  et  al.^  ^ 

All  calculations  of  the  optimum  controllers  in  the  presence 
of  noise  are  carried  out  in  the  domain.  This  means  the 

calculations  are  not  simple  in  form.  Since  the  method  involves 
the  inverse  jj__-  transformation,  it  is  difficult  to  write  the 
equations  as  part  of  a  digital  computer  program. 

Tou's  method  also  requires  a  specified  model  response 
to  be  incorporated  in  the  calculations.  No  such  model  can  be 
predicted  for  the  extended  deadbeat  controllers.  Because  of 
these  difficulties,  Tou's  method  is  not  considered  further. 

Another  method  considered  is  a  "straight  line" 
approximation.  This  means  that  the  deadbeat  response 
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obtained  using  the  extended  controller  is  given  by  an 

output  which  has  the  minimum  deviation,  at  the  sampling 

instants,  from  a  straight  line  from  c[0]  =  0  to  c[nT]  =  1.0. 

The  values  of  the  "straight  line"  output  at  the  sampling 

instants  are:  ccr  [iT]  =  -  ;  i  =  1,  2,  ...,  n-1.  Thus  the 

oil  n 

errors  for  the  straight  line"  at  the  sampling  instants  are: 

~  (1~  i/n)  ;  i  =  1,  2,  .  .  .  ,  n  — 1. 

For  the  extended  controllers,  the  deviations  from 
the  "straight  line"  approximation  are  given  by: 

i  . 

th 

Deviation  at  l  sampling  instant  e  DEV^gL  = 

(eiSL  ei}  ;  1=1/2,...,  n-1 

Therefore  the  sum  of  these  deviations  must  be  minimized  to 
yield  the  predicted  controller. 

This  method  is  considered  as  it  should  have  the 
smallest  change  in  the  equivalent  gain  of  the  controller 
from  one  period  to  the  next. 

The  proposed  criterion  was  applied  to  some  of  the 
plants.  It  was  found  to  be  a  useful  approximation  in  some 
cases,  but  the  results  in  general  were  not  satisfactory. 

Some  instances  when  this  is  a  useful  approximation  are 
presented  in  Chapter  7. 

The  failure  of  the  "straight  line"  approximation 
to  give  general  results  indicates  a  more  complicated  criterion 
is  needed.  A  more  general  criterion  is  derived  in  Section 


3-5. 
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3-5  DERIVATION  OF  THE  DESIGN  CRITERION 

In  the  system  shown  in  FIG.  1.1,  the  plant  parameters 
and  sampling  rate  are  assumed  to  be  fixed.  This  leaves  the 
parameters  of  the  extended  controller  as  the  only  variables. 
These  must  then  be  calculated  so  as  to  minimize  the  system 
noise  transfer. 

Since  the  plant  is  linear,  the  output  of  the  system 
at  any  sampling  instant  may  be  related  to  the  output  of  the 
controller  at  that  instant.  This  relationship  may  be  expressed 
as  follows: 

c[mT]  =  cG[mT]  BE2 [mT]  (3.1) 

where  c~[mT]  is  the  combined  equivalent  gain  of  the  plant 

th 

and  zero-order  hold  in  the  n\  sampling  period. 

This  expression  is  arrived  at  by  consideration  of  the 
discrete  state  equations  of  the  system.  In  these  equations, 
the  plant  is  replaced  by  a  matrix  of  values  which  vary  with 
time  but  are  constant  over  any  one  sampling  period. 

Using  (3.1),  the  mean  square  error  of  the  system,  as 
evaluated  at  the  sampling  instants,  becomes: 


Ic[mT]]2  =  [cQ [mT]  BE2 [mT ] ] 2  (3.2) 

Rewriting  (1.47)  as  a  difference  equation  yields  the 

x. 

controller  output  at  the  m  sampling  instant  as: 
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m 

BE 2 [mT]  =  Z  (h  )  BEl[mT  -  kT]  + 
k=0 

m 

Z  (e,  )  BE 2  [mT  -  kT]  (3.3) 

k-1 

where 

0  ^  m  ^  n  (3.4) 

Expanding  (3.3)  for  various  values  of  m  yields  the 
following  relationships. 

For  m  =  0 : 

BE 2  [ 0 ]  =  (hQ)  BE1 [ 0 ]  (3.5) 

For  m  =  1: 

BE2[T]  =  (  (hQ)  BEl  [T]  +  (h^)  BE1[0])  - 

((ex)  BE 2  [  0 ] )  (3.6) 

For  m  =  2 : 

BE2  [2T]  =  (  (hQ)  BEl  [ 2T ]  +  (h^)  BEl  [T]  + 

(h2)  BEl  [ 0 ] )  -  ((e1)  BE 2 [T]  - 

(e2)  BE 2  [  0 ]  ) 


(3.7) 


[Tm] saa 
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For  m  =  3: 

BE2 [ 3T]  =  (  (hQ )  BEl [ 3T]  +  (h±)  BE1[2T]  + 

(h2)  BEl [T]  +  (h3)  BEl [ 0 ] )  - 

(  (e1)  BE 2  [ 2T]  +  {e  )  BE2 [T]  + 

(e3)  BE 2  [ 0 ]  )  (3.8) 

This  expansion  can  be  similarly  extended  for  m  up  to  n. 

These  expressions  may  be  rewritten  entirely  in  terms  of  the 
input  BEl [mT]  in  the  following  fashion. 

For  m  =  0 : 

BE 2 [ 0 ]  =  (hQ )  BEl [ 0 ]  (3.5) 

For  m  =  1,  substituting  (3.5)  into  (3.6)  yields  in 
simplified  form: 

BE 2 [T]  =  (hx  -  e1hQ)  BE1[0]  +  (hQ)  BEl [T]  (3.9) 

For  m  =  2,  substituting  (3.5)  and  (3.9)  into  (3.7)  and 
simplifying  yields: 

BE 2  [2T]  =  (h2  -  ~  e2ho  +  elh0^  + 

(h±  -  ej^hQ)  BEl  [T]  +  (hQ)  BEl  [ 2T] 


(3.10) 


s  £  *  m  to*? 
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For  m  =  3,  substitute  (3.5),  (3.9)  and  (3.10)  into  (3.8) 

and  simplify  to  obtain: 


67 


BE 2 [3T] 


2e1e2hQ  -  e13hQ)  BEl[0]  +  (h2  - 


e2h0  +  el2ho}  BE1[T]  +  (hi  "  eiho} 


BE1 [ 2T]  +  (hQ )  BEl [ 3T] 


(3.11) 


As  shown  in  Table  3.3,  the  controller  inputs  BE1[0] 
to  BE1[3T]  are  acted  upon  by  the  same  set  of  controller 
coefficients.  This  set  of  coefficients  is  delayed  for  each 
input  BEl[mT]  by  mT  as  compared  to  BE1[0].  It  is  therefore 
possible  to  conclude  that  the  net  effect  of  the  controller 
over  all  n  periods  is  the  same  for  BE1[0]  as  for  any  other 
input . 

Rewrite  (3.5),  (3.9),  (3.10)  and  (3.11)  considering 

only  the  terms  with  respect  to  BE1[0]  as  follows: 


BE 2  [0] 


(hQ)  BE  [ 0 ] 


(3.5) 


BEl  [0] 


BE 2  [T] 


(hx  -  eih0)  BEl [ 0 ] 


(3.12) 


BEl  [0] 


BE 2 [2T] 


BEl  [0] 


BEl  [0] 


(3.13) 


tniBJ do  oS  Y3ilqrax8  bri£ 
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BE2 [3T] 


BE1 [0] 


(h3  e1h2  e2hi  e3ho  + 


el  hl  +  2ele2h0  el  h0 ^ 


BEl [ 0 ]  (3.14) 

If  the  net  controller  effect  on  any  BEl [mT]  is  constant, 
then  it  should  only  be  necessary  to  minimize  the  net  effect 
of  the  controller  on  BE1[0].  This  corresponds  to  minimizing 
the  effects  of  the  inputs  to  the  controller  at  all  sampling 
instants.  Since  the  error  signal  BEl [mT]  contains  the  noise 
signal,  it  is  implied  that  the  noise  transfer  through  the 
controller  is  minimized  by  this  calculation. 

This  means  that  the  "weighting  sequence"  BE2 [mT]  of 
the  extended  controller  is  minimized. 

Expressing  (3.1)  in  discrete  form  yields: 


c  [mT] 


BEl  [0] 


c  [mT]  BE 2 [mT] 

VJ 


BEl [0] 


(3.15) 


The  mean  square  error  due  to  BE1[0]  becomes: 


2 

c  [mT] 

= 

c  [mT]  BE2[mT] 

BEl [0] 

Cj 

• 

BEl [0] 

- 

(3.16a) 


Ignoring  the  constant,  (3.16a)  becomes: 


- 

2 

c[mT] 

= 

BE 2 [mT] 

BEl  [0] 

BEl  [0] 

(3.16b) 
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where  the  net  mean- square  error  of  the  controller  to  the 
input  BE1 [ 0 ]  is  defined  as: 


BE2 [mT] 


BE1 [0] 


f  1  ^ 
'BEl [ 0 ] ' 


n 

l 

m=0 


(BE 2  [mT] 


BE1  [0] 


(3.16c) 


Therefore  to  minimize  the  net  effect  of  the  controller  on 
BE1[0],  (3.16c)  must  be  minimized. 

Substituting  (3.5),  (3.12),  (3.13)  and  (3.14)  into 

(3.16c)  yields  the  following  simplified  form: 


BE 2  [mT] 


BE1 [0] 


(hQ) 2  +  (hx  -  e^hg) 2  + 


(h2  "  elhl  '  e2h0  +  el2h0)2  + 


(h3  -  e^h2  -  e2h1-  e3hQ  +  e^^  h1  + 


2ele2h0  el  h0)  +  *  *  *  + 


<BETT0T)2  ^[nT] 


) 


BEl [0] 


(3.17) 


This  sum  in  (3.17)  is  the  net  controller  mean- square  error 
due  to  BEl [ 0 ] .  Henceforth  this  term  will  be  referred  to  by 
the  abbreviation  NCE.  This  abbreviation  is  defined  as: 


s 
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NCE  =  NET  CONTROLLER  MEAN- SQUARE  ERROR 


TO  BEl [ 0 ] 


BE 2 [mT] 


BEl [0] 


(3.18) 


The  proposed  design  criterion,  which  when  minimized 
yields  the  predicted  extended  controllers,  is  obtained  by 
combining  (3.17)  and  (3.18)  into  the  following  form: 


NCE  =  (hQ)2  +  (h±  -  exh0)2  +  (h2  -  e1h1  -  e^  + 


el  h0^  +  ^h3  elh2  e2hl  e3h0  +  el  hl  + 


2ele2h0  el  h0)  +  * 


+  /  1. 
BEl  [  0]  ; 


(BE 2 [nT] 


) 


BEl  [0] 


(3.19) 


This  criterion  is  defined  entirely  in  terms  of  the 
parameters  of  the  controller.  These  parameters,  in  turn,  are 
functions  of  the  plant  parameters  and  the  sampling  rate  of  the 
system  only. 

Therefore  even  if  a  particular  plant  has  a  peak  or 
dip  at  some  frequency,  the  extended  controller  predicted  by 
NCE  should  still  give  the  minimum  noise  transfer  for  the 
bandlimited  white  noise.  This  conclusion  is  based  on  the 
fact  that  NCE  may  be  written  entirely  as  a  function  of  the 
plant  parameters  and  the  sampling  rate.  Therefore,  any  dips 
or  peaks  in  the  frequency  response  of  the  plant  will  be 
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considered  in  the  calculation  of  the  predicted  extended 
controller . 

3-6  COMPARISON  OF  CRITERION  RESULTS  AND 
HYBRID  RESULTS 

For  the  five  systems  considered  at  the  start  of 
this  chapter,  the  criterion  of  section  3-5  is  used  to  find 
the  values  of  e^  which  yield  the  minimum  values  of  NCE. 

These  results  are  compared  in  Table  3.4  to  the  hybrid  results 
listed  in  Table  3.2. 

From  the  results  in  Table  3.4,  it  is  apparent  that 
acceptable  agreement  exists  between  the  hybrid  and  NCE  values 
of  e^.  In  all  cases,  the  value  of  e^  predicted  by  NCE  lies 
in  the  middle  of  the  range  of  e^  values  from  the  hybrid. 

These  results  indicate  that  this  proposed  criterion 
is  worth  additional  consideration.  To  further  test  this 
criterion,  more  simulation  results  for  the  predicted 

controllers  are  required. 

In  Chapter  4,  an  experimental  simulation  using  a 

digital  computer  is  developed  to  supply  these  additional 


results . 


PLANT  (SECONDS) 
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TABLE  3.4  COMPARISON  OF  HYBRID  AND  NCE  RESULTS 
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CHAPTER  IV 


DIGITAL  COMPUTER  SIMULATION 

In  this  chapter,  the  digital  computer  simulation 
used  to  test  the  controllers  calculated  by  NCE  is  outlined. 
The  hybrid  simulation  results,  with  which  the  predictions 
by  the  design  criterion  are  compared  in  Section  3-6,  are 
for  only  two  types  of  plants.  The  digital  simulation  is 
used  to  test  the  four  types  of  plants  listed  in  Chapter  I. 
The  results  of  these  simulations  are  given  in  Chapters  5 
to  7 . 

The  digital  simulation  is  used  instead  of  the  hybrid 
simulation  for  the  following  reasons: 

1.  The  controllers  can  be  tested  more  rapidly. 

2.  A  larger  combination  of  plants  and  controller 
extensions  may  be  considered. 

3.  Repeatable  results  are  obtainable  with  precisely 
defined  noise  sources. 

4.  The  simulation  is  free  of  the  component  drift 
present  in  the  hybrid  simulation. 

5.  Accuracy  of  test  results  up  to  16  significant 
figures  is  possible  if  desired.  The  actual 
calculation  accuracy  is  5  significant  figures. 

In  the  first  section,  the  method  of  simulation  is 
outlined.  Then  the  test  procedure  followed  is  given.  The 
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noise  sources  used  in  the  test  program  are  discussed  in  the 
final  section  of  this  chapter. 


4-1  METHOD  OF  DIGITAL  SIMULATION 

Fortran  IV  is  the  programming  language  used  on  an 
IBM  360/67  computer.  The  systems  are  simulated  by  use  of 
the  discrete  state  equations  derived  in  Chapter  1. 

Various  noise  sources  are  simulated  by  the  random 
number  generation  programs  available  in  the  IBM  Scientific 
Subroutine  Package.  Since  the  mean-square  error  of  the 
system  is  calculated  only  at  the  sampling  instants,  the  input 
noise  need  only  be  defined  at  the  sampling  instants  to  yield 
the  system  noise  transfer  value. 

The  mean-square  sampled  error  may  be  calculated 
theoretically  for  systems  with  a  noise  contaminated  input. 

To  perform  this  calculation,  it  is  necessary  to  know  the 
auto-spectral  and  cross-spectral  densities  of  the  unit  step 
and  noise  inputs.  It  is  also  necessary  to  take  the  inverse 

-  transform  of  the  system  transfer  function.  The  variable 
s  is  then  replaced  in  the  resultant  function  by  -s  and  the 
transform  is  then  taken.  Then  a  contour  integral  is 
evaluated. 

This  calculation  is  expressed  in  equation  form  by 


Tou1  as: 
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(nT) 


(4.1) 


where  <J>  (z)  =  M(z  1)  M(z)  <|>  (z) 

ee  Ynn 


(4.2) 


and  where  M(z)  =  C(z)/R(z) 

In  this  calculation,  the  following  assumptions  are  made: 


1.  There  is  no  correlation  between  the  noise  and 
the  unit  step  input. 

2.  The  model  system  response  is  set  equal  to  the 
system  response  obtained  using  the  extended 
controller  being  tested.  This  removes  the  terms 
from  (4.2)  which  contain  the  auto-spectral  density 
of  the  unit  step  d>  (z)  . 

3.  4>nn(z)  is  t*le  auto-spectral  density  of  the  input 

noise.  For  white  noise,  (j>  (z)  is  a  constant. 

These  calculations  are  very  difficult  to  program  on 


the  digital  computer  due  to  the  transforms  and  contour 
integration  required.  Therefore  the  alternate  approach  of 
using  random  number  sets  to  simulate  the  white  noise  input 
is  used. 


This  method  yields  only  an  approximation  to  white 


noise  as  the  experiment  is  time  limited.  The  validity  of 
any  results  obtained  by  this  method  may  be  evaluated  if  the 
accuracy  of  the  approximation  is  known.  These  random  number 
approximations  are  discussed  in  Section  4-3. 


■  :  a  ; 
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The  Fortran  program  consists  of  the  following 
sections : 

1.  A  group  of  algorithims  to  calculate  the  various 
controllers . 

2.  A  system  simulation  to  verify  that  deadbeat  to 
a  unit  step  is  satisfied. 

3.  A  system  simulation  to  test  the  system  in  the 
presence  of  noise. 

4.  An  algorithim  to  calculate  the  mean-square  error 
for  the  minimum  time  controllers  and  the  normalized 
mean-square  errors  of  the  extended  controllers. 

The  block  diagram  of  the  simulation  given  by  FIG.  4.1 
incorporates  these  sections. 

The  diagram  of  the  computer  simulation  shown  in 
FIG.  4.1  is  greatly  simplified.  All  details  are  omitted. 

The  program  for  the  sections  which  verify  deadbeat  and 
calculate  the  noise  transfers  is  given  in  Appendix  B. 

The  algorithims  for  the  calculation  of  the  various 
predicted  extended  controllers  are  given  in  the  appropriate 
chapters . 

4-2  DIGITAL  COMPUTER  SIMULATION  TEST  PROCEDURE 

As  shown  in  FIG.  4.1,  there  are  three  types  of 
controllers  to  be  considered  for  a  given  system.  The 
computer  simulation  package  can  be  run  in  two  forms.  These 


are : 
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FIGURE  4.1  SIMPLIFIED  BLOCK  DIAGRAM  OF  DIGITAL  COMPUTER  SIMULATION 
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a.  Branches  (1)  and  (2)  compared. 

b.  Branches  (1)  and  (3)  compared. 

Both  forms  begin  by  calculating  the  minimum  time  controller 
and  its  corresponding  noise  transfer.  Then  in  form  (a),  the 
predicted  controller  is  calculated  and  its  noise  transfer  is 
normalized. 

In  form  (b) ,  the  error  terms  of  the  extended  controller 
are  iterated  as  in  Section  2-4.  This  yields  a  series  of 
controllers  which  are  then  tested  in  the  presence  of  noise 
to  find  the  minimum  normalized  noise  transfer  controller. 

This  controller  is  called  the  "grid  searched"  controller. 

This  "grid  search"  is  similar  to  the  test  procedure  used  in 
the  hybrid  computer  simulation  described  in  Chapters  2  and  3. 

"Grid  searched"  controllers  are  calculated  for  various 
combinations  of  plants,  sampling  rates  and  controller 
extensions.  These  controllers  are  then  compared  with  the 
appropriate  predicted  controller.  From  these  comparisons, 
the  accuracy  of  the  prediction  by  the  design  criterion  is 
obtained.  Results  of  these  tests  are  presented  in  Chapters 
5 ,  6  and  7 . 

4-3  NOISE  SOURCES  SIMULATED  BY  THE  RANDOM 
NUMBER  GENERATORS 

It  was  stated  earlier  that  the  input  noise  for  the 
system  is  assumed  to  be  white.  However  it  is  not  possible 
in  a  time  limited  experiment  to  obtain  wide-band  white  noise. 

It  is  possible  to  generate  noise  which  is  white  over  some 
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limited  range  of  frequencies.  This  is  called  bandlimited 
white  noise. 

After  sampling,  the  noise  bandwidth  is  reduced  to  0 

to  co  /2  where  co  is  the  sampling  frequency.  If  the  noise  in 
s  s 

this  frequency  range  is  white,  the  distribution  of  frequencies 

is  uniform  from  0  to  cog/2 .  This  occurs  if  the  input  noise 

before  sampling  is  white  and  tightly  low-passed  from  0  to 

co  /2.  Then  the  output  of  the  sampler  is  pure  bandlimited 
s 

white  noise.  In  this  case,  the  bandlimited  noise  referred  to 
is  actually  low-pass  noise  in  all  cases. 

It  is  useful  to  calculate  the  Fourier  coefficients  of 
the  sampled  noise  sources  simulated  by  the  random  number  sets 
used  in  the  digital  simulation.  The  distribution  of  these 
coefficients  may  be  compared  to  the  ideal  bandlimited  white 
noise  frequency  distribution  obtained  after  sampling. 

Using  a  Fortran  subroutine  called  PS301A,  the  real 
and  imaginary  Fourier  coefficients  are  calculated  for  the 

noise  sources  by  the  Fast  Fourier  method.  The  frequencies 

co 

for  which  these  coefficients  are  calculated  are  0,  (75^-), 

Ws  °s  L-l 

2  (-7r=r)  ,  3(757-)  ,  .  .  .  ,  (-75: — )  co  where  L  is  the  total  number 

ZLi  ZLi  Z  s 

of  points  in  the  random  number  set.  At  each  point,  the 
magnitude  of  the  Fourier  coefficient  is  calculated  in  decibels 
from  the  real  and  imaginary  values. 

Once  these  magnitudes  are  known,  a  comparison  may  be 
made.  The  mean  of  these  magnitudes  for  each  noise  source  is 
calculated.  Then  the  variance  and  deviation  of  the  magnitudes 
around  the  mean  are  evaluated  for  each  source. 


. 
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Since  low-pass  white  noise  has  a  uniform  distribution 

over  0  to  oj  /2 ,  its  variance  and  deviation  around  the  mean 
s 

are  zero.  Therefore  the  lower  the  values  of  the  variance 
and  deviation  for  each  noise  source,  the  better  the  approxi¬ 
mation  . 

The  Scientific  Subroutine  Package  of  Fortran  IV 
contains  two  random  number  generation  programs.  These  are 
called : 

1 .  RANDU 

2.  GAUSS 

The  subroutine  RANDU  generates  uniformly  distributed 
random  numbers  between  0  and  1.0.  In  order  to  make  this 
sequence  useful,  0.5  is  subtracted  from  all  values  when 
calculated.  This  results  in  a  uniformly  distributed  set  of 
numbers  with  amplitudes  between  -0.5  and  +  0.5. 

The  second  subroutine  GAUSS  generates  a  set  of 
random  numbers  with  a  specified  variance  around  a  specified 
mean.  The  amplitude  distribution  of  these  numbers  is 
Gaussian.  For  all  Gaussian  sets  used,  the  mean  is  0,  but 
the  variances  differ. 

By  varying  initializing  coefficients  in  the 

subroutines,  different  sets  of  random  numbers  are  available. 

The  length  of  the  set  of  random  numbers  is  variable  up  to 
29 

10  numbers  without  repetition  of  values.  For  normal  testing, 
the  sets  generated  contain  256  or  1024  numbers.  The  reasons 
for  these  lengths  are  given  in  Chapter  5. 
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Five  different  noise  sources  are  used.  These  are 

labeled : 

1.  GAUSS 2 A 

2.  GAUSS 2 

3.  GAUSS 3 

which  are  generated  by  the  GAUSS  subroutine  and 

4 .  RANDUA 

5 .  RANDU 

which  are  generated  by  the  RANDU  subroutine. 


Table  4.1 

lists  the 

Fourier  statistics 

calculated  for  each 

noise  source  by  the 

subroutine  PS301A. 

NOISE 

LENGTH 

OF 

FOURIER  COEFFICIENT 

SOURCE 

NUMBER 

SET 

STATISTICS 

MEAN 

DEVIATION  ABOUT 

(decibels ) 

MEAN  (decibels 

GAUSS 2A 

1024 

7.6053 

5.4287 

GAUSS 2 

256 

1.4664 

5.2050 

GAUSS 3 

256 

1.7052 

5.3203 

RANDUA 

1024 

16.746 

5.7010 

RANDU 

256 

10.618 

5.6861 

TABLE  4.1 

FOURIER 

COEFFICIENT  STATISTICS 

FOR  RANDOM 

NUMBER  NOISE  SOURCES 
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In  Table  4.1,  only  the  deviation  is  given  as  the 
variance  is  simply  the  square  of  the  deviation.  From  this 
table,  GAUSS2  appears  to  be  the  best  approximation  as  it  has 
the  lowest  deviation.  In  all  cases,  the  deviation  is 
substantially  greater  than  zero  however. 


c  :  n  -  Icidu* 


CHAPTER  V 


CONTROLLERS  EXTENDED  ONE  TERM 

The  results  of  the  digital  simulation  tests  of  the 
controllers  with  one  extra  term  are  presented  in  this 
chapter.  To  simplify  the  presentation,  a  controller 
extended  by  one  term  is  referred  to  as  a  PIC.  This  stands 
for  Plus  One  Term  Controller. 

The  calculation  algorithim  for  the  PIC  is  derived 
from  NCE  in  the  first  section.  In  the  following  sections, 
the  test  results  for  the  four  plants  listed  in  Chapter  1 
are  presented. 

The  results  for  these  plants  consist  of  two  parts. 

1.  The  controllers  predicted  by  NCE  are  presented 
for  various  sets  of  system  parameters. 

2.  The  predicted  controllers  are  compared  to  the 
"grid  searched"  controllers  using  the  various 
noise  sources. 

The  responses  of  the  minimum  time  deadbeat  controllers 
and  the  predicted  PlCs  to  ramp  and  acceleration  inputs  are 
compared  in  the  second  last  section.  In  the  last  section, 
some  conclusions  obtained  from  the  test  results  are  given. 
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The  four  plants  defined  by  (1.4)  to  (1.7)  were  tested 
using  digital  simulation.  In  order  to  test  these  plants, 
the  PlCs  must  be  calculated.  The  method  presented  in 
Section  2-4  is  used  as  the  basis  of  the  algorithim. 

Consider  the  plant  G2(s).  In  Table  2.2(b)  the 
coefficients  of  the  D(z)  are  expressed  in  terms  of  e^. 

Since  NCE  is  derived  in  Chapter  3  in  terms  of  the 
coefficients  of  the  D(z),  it  is  possible  to  rewrite  NCE  as 
a  function  of  e^. 

The  calculation  algorithim  yields  the  predicted 
value  of  e-^  using  NCE.  This  allows  the  predicted  PIC  to  be 
calculated.  Although  G 2 (s)  is  used  as  an  example  during 
this  derivation,  the  algorithim  derived  is  general. 

The  digital  controller  coefficients  listed  in  Table 
2.2(b)  may  be  expressed  as  follows: 


h 

h 

h 

h 

e 


0 


1 


2 


3 


2 


alel  +  a2 


=  a3ei  +  a4 


a5el  +  a6 


=  a7el  +  a8 


=  b  e±  +  b2 


(5.1) 

(5.2) 

(5.3) 

(5.4) 


(5.5) 
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e3  =  b3e1  +  b4  (5.6) 

where  the  coefficients  to  ag  and  to  for  G 2  (s)  are 
defined  in  Table  5.1(a). 

The  criterion  proposed  in  Chapter  3  for  extended 
controllers  with  n  ^  3  is 


NCE  =  hQ2  +  (h1  -  hoe]_)2  +  (h2  “  hlel  "  h0e2  + 

h0el2)2+  (h3  •  h2el  -  hle2  -  h0S3  +  hlel2  ' 
h0Sl  +  2h0Sle2*  ’  (5.7) 

(5.7)  may  be  rewritten  in  the  following  form: 

n 

NCE  =  Z  NCE.  (5.8) 

i  =  0  1 


where 


NCE0  = 

<v2 

NCE1  = 

- 

,  x  2 
hoei 

nce2  = 

(h2 

hlel  -  hoe 

NCEg  = 

(h3 

h2Sl  _  hlS 

.  3 
h0ei 

+  2hoeie2) 

2,2 


l0cl 


0  3  11 


(5.9) 

(5.10) 

(5.11) 


(5.12) 
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For  G2(s),  the  PIC  has  n  =  3.  By  substituting  (5.1) 


to  (5.6)  into  (5, 

.9)  to  (5.12),  NCE  may  be  written  as  a 

function  of  e^. 

Define  NCE  for  a  PIC  as  NCE^  which  is  a 

function  of  e^. 

Using  this  definition  yields  (5.9)  as: 

NCEio  = 

2  2  2 

al  el  +  ^ala2el  +  a2  (5.13) 

(5.10)  as : 


nce1jl  = 

2  4  o  3  ,  2  0  x 

al  el  2alallel  +  ^all  2ala4 

e^2  +  2a4a11e1  +  a42  (5.14) 

(5.11)  as: 


NCE12  “ 

al  el  +  2ala21el  +  (a21  +  2al  a22 ^ 

el  +  2 ^ala23  +  a21a22)el  +  ^a22  + 

2a21a23^el  +  2a22a23Sl  +  a23  (5.15) 

(5.12)  as: 


nce13  = 

2  3  7  2  6 

al  el  “  2ala31el  +  (a31  “  2ala32)el  + 

2 (a31a32  ”  ala33)el  +  ^a32  +  2a31a33 

2ala34)el  +  2(a31a34  +  a32a33)el  +  (a33  + 

2a32a34)el  +  2a33a34el  +  a34  (5.16) 

(5.16) 
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where 


all  = 

a3  "  a2 

a21 

a2  -  a3  -  a1b1 

a22 

a5  ~  a4  alb2  ~  a2bl 

a2  3 

3.c  —  a~b~ 

6  2  2 

a31 

a3  -  a2  +  2a1b1 

a32 

a4  -  a5  -  a3b1  -  a1b3  +  2  (a^  +  a^) 

a33 

a_  -  a,  +  2anbn  -  a~bn  -  a„b,  -  a,b.  -  anb_ 

7  6  22  32  41  14  23 

a34 

a8  ~  a4b2  “  a2b4  (5.17) 

Substituting  (5. 

.13)  to  (5.16)  into  (5.8)  yields  the  following 

simplified  form: 


nce1  = 

2n+2  .  n 

Z  C.e  1  =  Z  (NCE  . )  (5.18) 

i=0  1  1  i=0  11 

where  C.  is  the 

l 

coefficient  of  e^1  generated  by  summing  the 

coefficients  of  e^1  in  (5.13)  to  (5.16).  Thus  for  G2(s),  the 


criterion  is  an 

th 

8  order  linear  equation  in  e^  as  shown  in 

(5.19)  . 


NCE1  = 

8 

Z  C.e  1  (5.19) 

i=0 

. 

neioiiidoo 
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This  expression  must  be  minimized  in  order  to  obtain  the 
predicted  value  of  e^.  Because  of  the  deadbeat  restriction, 
e^  for  any  valid  controller  must  satisfy  the  following 
conditions: 

e ^  is  real  (5.20a) 

0  -  e±  -  1.0  (5.20b) 

One  method  for  finding  the  minimum  of  an  equation  is 
differentiation.  Differentiating  (5.19)  with  respect  to  e^ 
yields : 


(NCE1) 


8 

L 

i=l 


i 


C.e 

l 


i-1 

1 


(5.21) 


Set  (5.21)  to  zero: 


(NCE1) 


0 


8 

l 

i=l 


i 


C.e 

l 


i-1 

1 


(5.22) 


The  roots  of  this  equation  yield  the  values  of  e^  which 
correspond  either  to  the  local  maxima  or  minima  of  the 
expression  NCE^ . 

A  polynomial  root  solving  subroutine  called  POLRT 

is  available  in  the  Fortran  Scientific  Subroutine  Package. 

This  solves  for  the  real  and  imaginary  parts  of  all  the 

th 

roots  of  a  polynomial  of  up  to  36  order. 

Any  roots  not  satisfying  the  restrictions  of  (5.20) 
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may  be  ignored.  It  was  found  that  for  every  combination  of 
a,  3  and  T  checked,  only  one  root  of  (5.21)  would  satisfy 
(5.20).  For  these  cases,  the  calculated  value  of  e^ 
corresponded  to  the  only  minimum  of  NCE^. 

Similar  calculations  were  performed  for  G^(s), 

Gg(s)  and  G^(s).  For  all  combinations  of  parameters 
attempted,  only  one  value  of  e^  was  found  in  each  case  which 
satisfied  (5.20).  In  each  case  this  value  of  e^  corresponded 
to  the  only  minimum  of  NCE^. 

The  coefficients  a1  to  ag  and  b-^  to  b^  are  listed  in 
Table  5.1(a)  and  5.1(b)  for  G-^(s)  to  G^(s).  These  coeffi¬ 
cients  are  defined  in  terms  of  the  appropriate  constants  from 
Chapter  1. 

From  these  observations,  a  general  algorithim  for  the 
calculation  of  the  predicted  PIC  is  as  follows: 


1.  Determine  n  which  is  the  order  of  the  PIC  to  be 
calculated. 

2.  Calculate  NCE^  as  shown  in  (5.18). 

2n+2 

NCE  =  I  C.e  1  (5.18) 

1  i=0  1 

3.  Differentiate  (5.18)  with  respect  to  e^  and 
set  the  result  to  zero. 


(NCEX) 


2n+2  . 

=  I  i  C1e11  =  0 

i=l 


(5.23) 


Solve  for  the  roots  of  (5.23)  and  choose  as 
e^  the  root  which  satisfies  (5.20). 


4. 


,  *  . 
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5.  Substitute  this  value  of  e^  into  (5.1)  to  (5.6) 

to  calculate  the  predicted  minimum  noise  transfer 
PIC. 

5-2  RESULTS  FOR  THE  PICs  of  G^s)  =  1/s  (s+«) 


The  greatest  amount  of  work  was  done  for  this  plant 
as  it  was  the  first  and  simplest  to  be  tested.  Also  any 
conclusions  obtained  from  it  may  be  applicable  to  later 
plants  tested.  The  test  procedure  followed  is  outlined  in 
Chapter  4. 

The  sampling  rate  T  and  the  pole  a  are  the  only 
system  parameters  considered  in  NCE^.  Therefore  T  was  fixed 
and  successive  values  of  a  were  used  in  testing  the  plants. 

The  first  set  of  data  collected  for  G-^(s)  is 
obtained  using  the  GAUSS2A  and  RANDUA  noise  sources.  Both 
are  1024  sampling  periods  long.  FIG.  5.1  shows  the  values 
of  e^  predicted  by  NCE^  plotted  versus  «  with  T  =  1  second. 

It  is  only  necessary  to  present  e^  as  all  the  other 
coefficients  of  the  D(z)  can  be  calculated  using  the 
relationships  given  in  Table  5.1. 

FIG.  5.2  presents  the  normalized  noise  transfers  of 
the  systems  using  the  predicted  PICs  when  the  noise  sources 
are  GAUSS2A  and  RANDUA.  The  noise  transfer  is  shown  in  all 
following  graphs  as  a  fraction.  This  means  the  noise  transfer 
of  the  system  containing  the  PIC  is  normalized  to  the  noise 
transfer  of  the  system  containing  the  appropriate  minimum 
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time  deadbeat  controller. 

All  calculations  following  are  carried  out  for 
T  =  1  second  as  the  predicted  values  of  e^  are  constant  for 
constant  values  of  the  normalized  constant  aT.  For 
example  if  T  =  10  seconds  and  =  0.1,  the  e^  predicted  is  the 
same  as  for  «  =  10  and  T  =  0.1  second. 

A  grid  search  was  then  performed  with  both  GAUSS2A 
and  RANDUA.  The  number  of  values  of  «  considered  is  less 
than  for  the  predicted  PICs  as  the  computer  time  to  perform 
the  grid  search  is  much  longer.  Also  a  good  indication  of 
the  agreement  should  be  obtained  by  taking  a  reasonable 
cross-section  of  values.  In  order  to  further  reduce  the 
computer  time  required,  the  value  of  e^  was  only  calculated 
to  two  decimal  places.  FIG.  5.3  presents  the  comparison 
between  the  predicted  values  of  e^  and  the  values  of  e^ 
found  by  the  grid  searches. 

The  percentage  differences  between  the  predicted  and 
grid  searched  PICs  are  now  considered.  All  the  differences 
between  the  D(z)  error  terms  are  calculated  as  a  percentage 
of  the  predicted  value.  The  normalized  noise  transfer 
differences  are  expressed  as  a  percentage  of  the  predicted 
noise  transfer. 

The  noise  transfer  difference  using  GAUSS2A  is 
-  1.5%  at  ttT  =  0.1  and  zero  at  aT  =  96.0.  Using  RANDUA, 
the  difference  is  -  3.3%  at  «T  =  0.1  and  zero  at  aT  =  96.0. 

The  comparable  differences  in  e^  are  -  4.4%  and  +  2.6%  from 
the  GAUSS 2 A  grid  search.  Using  RANDUA,  the  differences  are 
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FIGURE  5.1  PREDICTED  e,  FOR  PlCs  OF  G,  (s)  WITH  T=1  SEC 


NORMALIZED  NOISE  TRANSFER 
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0.1  1.0  10.0  100.0 

a 

FIGURE  5.2  NORMALIZED  NOISE  TRANSFERS  USING  PREDICTED 


PlCs  OF  G^s)  WITH  T=1  SEC.  USING  GAUSS2 
AND  RANDUA 

FIGURE  5.3  PREDICTED  AND  GRID  SEARCHED  USING 

GAUSS 2 A  AND  RANDUA,  FOR  PlCs  OF  G^s) 
WITH  T=1  SEC. 
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-  5.5%  at  «T  =  0.1  and  -  0.9%  at  aT  =  96.0.  This  agreement 
is  considered  acceptable. 

The  results  obtained  to  this  point  were  generated  by 
test  runs  of  1024  sampling  periods.  This  is  lengthy  in  terms 
of  computer  time.  Shorter  runs  were  attempted  to  determine 
if  the  results  are  comparable.  The  length  was  reduced  to 
256  sampling  periods.  The  new  noise  sources  are  GAUSS2, 

GAUSS 3  and  RANDU. 

FIG.  5.4  is  a  plot  of  the  mean  of  the  normalized 
noise  transfers  calculated  for  all  the  noise  sources  at  each 
value  of  «.  The  95%  confidence  limits  are  shown  with  each 
mean . 

Using  these  noise  sources,  grid  searches  were  carried 
out.  The  95%  confidence  limits  of  the  mean  values  of  e^ 
calculated  from  these  grid  searches  are  compared  to  the 
predicted  values  of  e^  in  FIG.  5.5. 

Table  5.2  shows  the  differences  for  the  normalized 
noise  transfers  calculated  using  the  different  noise  sources. 
The  differences  between  the  calculated  and  grid  searched 
values  of  e^  are  also  shown. 
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a 


FIGURE  5.4  95%  CONFIDENCE  LIMITS  ON  MEAN  NORMALIZED 

NOISE  TRANSFER  USING  PlCs  OF  G, (s) 

WITH  T=1  SEC.  1 

FIGURE  5.5  95%  CONFIDENCE  LIMITS  ON  MEAN  VALUES  OF 

GRID  SEARCHED  e±  COMPARED  TO  PREDICTED  e± 

FOR  PlCs  OF  G1(s)  WITH  T=1  SEC. 
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ei 

Differences 

Normalized  Noise 

Transfer  Differences 

Noise  Source 

«  =  0.1 

«  =  96 

.0 

a  =  .  1 

=  96.0 

GAUSS 2 A 

-4.4% 

+  2.6% 

-1.5% 

NIL 

RANDUA 

-5.5% 

-0.9% 

-3.3% 

NIL 

GAUSS 2 

-1.6% 

+  6.4% 

-0.2% 

-0.3% 

GAUSS 3 

-5.4% 

+  2.7% 

-3.3% 

-0.1% 

RANDU 

-6.8% 

+  0.9% 

-3.7% 

NIL 

TABLE  5.2 

COMPARISON  OF 

PREDICTED 

AND 

GRID  SEARCHED 

PlCs  FOR  G1(s)  WITH  T  =  1  SECOND 


From  Table  5.2,  it  is  evident  that  the  noise  sources 
all  give  comparable  results.  The  maximum  noise  transfer 
difference  is  3.7%  and  the  maximum  e^  difference  is  6.8%. 
This  indicates  that  shorter  test  runs  are  similar  in  results 
to  the  longer  runs.  Because  of  this  agreement,  the  shorter 
runs  are  used  for  most  of  the  plants  tested  later. 

5-3  RESULTS  FOR  THE  PlCs  of  G2(s)  =  1/s  (s+«)  (s+3) 

The  PlCs  of  G2(s)  have  n  =  3  and  the  minimum  time 
deadbeat  has  n  =  2. 

For  this  plant,  the  test  procedure  consisted  of 
fixing  various  values  of  6  and  T  and  then  varying  «  over  the 
range  0.1  to  96.0.  The  values  of  e^  as  calculated  by  NCE^ 
are  plotted  for  the  following  sets  of  parameters: 
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1.  3  =  0.1,  T  =  1  second  in  FIG.  5.6 

2.  3  =  1.0,  T  =  1  second  in  FIG.  5.6 

3.  3  =  10.0,  T  =  1  second  in  FIG.  5.7 

4.  3  =  100.0,  T  =  1  second  in  FIG.  5.7 

5.  3  =  1.0,  T  =  10  seconds  in  FIG.  5.7 

6.  3  =  1.0,  T  =  100  seconds  in  FIG.  5.7 

The  normalized  noise  transfers  for  the  systems 
containing  the  predicted  PICs  were  then  calculated  using 
GAUSS 2A.  This  noise  source  was  used  when  calculating  the 
noise  transfer  of  the  predicted  PICs.  The  shorter  test  runs 
were  used  in  the  grid  searches.  The  plants  for  which  the 
predicted  normalized  noise  transfer  is  plotted  are: 

1.  3  =  0.1,  T  =  1  second  in  FIG.  5.8 

2.  3  =  1.0,  T  =  1  second  in  FIG.  5.8 

3.  3  =  10.0,  T  =  1  second  in  FIG.  5.8 

4.  3  =  100.0,  T  =  1  second  in  FIG.  5.8 

As  in  Section  5.2,  a  grid  search  was  carried  out  at 
selected  values  of  the  plant  coefficients  to  compare  the 
predicted  and  grid  searched  values  of  e..  The  noise  source 
used  was  GAUSS2.  The  plant  used  had  3  =  1.0  and  T  =  1 
second  with  «  =  .1,  1.5,  12.0  and  96.0.  FIG.  5.9  shows  the 
values  of  e^  compared. 

The  comparative  results  are  always  within  0.5%  of 
each  other  for  the  noise  transfer  calculations.  Also  the 
agreement  between  the  predicted  and  grid  searched  controllers 
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FIGURE  5.6  PREDICTED  e1  FOR  PlCs  OF  G2(s)  WITH  T=1  SEC. 

FIGURE  5.7  PREDICTED  e,  FOR  PlCs  OF  G  (s)  WITH  T=1  & 

10  SEC. 
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a 

FIGURE  5.8  NORMALIZED  NOISE  TRANSFERS  USING  PREDICTED 

PICs  OF  G2(s)  WITH  T=1  SEC.  USING  GAUSS2A 


PREDICTED  e±  AND  GRID  SEARCHED  USING  GAUSS2 , 

FOR  PICs  OF  G2(s)  WITH  3=1  AND  T=1  SEC. 


FIGURE  5.9 
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is  best  at  a  =  0.1  for  both  the  noise  and  e^.  The  percentage 
differences  for  e^  are  -0.2%  at  «  =  0.1  and  -0.3%  at  «  =  96.0. 
These  results  again  are  considered  satisfactory. 

5-4  RESULTS  FOR  THE  PlCs  OF  G3 (s)  =  (s+3)/s(s+«) 

This  plant  was  tested  using  the  same  procedures  as 
in  Sections  5-2  and  5-3.  The  predicted  values  of  e^  were 
calculated  for  fixed  3  and  T  with  «  varied.  Predicted 
values  of  e^  for  the  PlCs  are  given  for  3  =  2.0,  5.0  and 
50.0  with  T  =  1  second  in  FIG.  5.10.  The  normalized  noise 
transfers  for  the  systems  containing  these  controllers  are 
calculated  using  GAUSS2  and  are  shown  in  FIG.  5.11. 

Following  the  same  procedure  used  previously,  a 
grid  search  using  GAUSS2  was  done  at  various  values  of  «  for 
3  =  5.0  and  T  =  1  second.  The  results  from  this  are  shown 
in  FIG.  5.12. 

The  differences  between  the  e^  values  are  -2.7%  at 
<*  =  0.1  and  +7.1%  at  a  =  96.0.  For  the  noise  transfer, 
these  differences  are  +0.3%  and  -0.4%  respectively.  Although 
the  e^  values  differ  by  quite  large  margins,  the  noise 
transfers  are  very  close.  This  indicates  that  the  change  in 
the  noise  transfer  is  very  slow  in  the  region  surrounding  the 
predicted  values  of  e^^  when  using  GAUSS 2 .  This  experimental 
agreement  is  considered  satisfactory. 

For  some  combinations  of  a,  3  and  T,  these  plants 
exhibit  an  overshoot  in  the  deadbeat  response  to  a  unit  step. 
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FIGURE  5.10  PREDICTED  FOR  PICs  OF  G3(s)  WITH  T=1  SEC. 
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a 

FIGURE  5.11  NORMALIZED  NOISE  TRANSFERS  USING  PREDICTED 

PICs  OF  G3(s)  WITH  T=1  SEC.  USING  GAUSS2 

FIGURE  5.12  PREDICTED  e  AND  GRID  SEARCHED  USING 

GAUSS 2 ,  FOR  PICs  OF  G3(s)  WITH  3=5  AND  T=1  SEC. 
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FIGURE  5.13  MAXIMUM  OVERSHOOTS  TO  A  STEP  INPUT  FOR  PREDICTED  PlCs  AND 

MINIMUM  TIME  DEADBEAT  CONTROLLERS  OF  G~(s)  WITH  T=1  SEC. 
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This  is  illustrated  in  FIG.  1.2(b).  The  minimum  time 
controller  and  the  PIC  both  exhibit  this  peak.  However  the 
PIC  reduces  the  height  of  this  peak  and  therefore  is  an 
improvement. 

As  an  illustration  of  this  improvement,  the  deadbeat 
responses  for  the  PICs  and  minimum  controllers  were  calculated 
for  3  =  2.0  and  5.0  with  T  =  1  second.  The  magnitudes  of 
these  peaks  are  plotted  versus  «  in  FIG.  5.13.  From  this 
graph,  the  peaks  for  the  PICs  are  only  half  as  high  as  the 
peaks  for  the  minimum  time  controllers . 

5-5  RESULTS  FOR  THE  PICs  OF  G4(s)  =  1/s  (s2+ 

2  ^o)  s  +  (i)  2) 

n_ n 

The  fourth  plant  considered  is  G4(s)  which  may  have 

complex  poles  in  the  s  -  plane  for  particular  plant 

parameters.  The  values  of  C  normally  are  between  0  and  1. 

To  test  this  plant  the  values  of  C  and  T  were  fixed 

and  a)  was  varied  from  0.1  to  96.0.  In  FIG.  5.14,  the  values 
n 

of  e^  predicted  from  NCE^  are  shown  for  T  =  1  second  and 

C  =  0.3,  0.5  and  0.7.  The  normalized  noise  transfers  of  the 

systems  containing  these  plants  when  subjected  to  GAUSS2  are 

shown  in  FIG.  5.15.  FIG.  5.15  shows  the  predicted  values  of 

e^  for  C  =  0*3  and  0.5  with  T  =  5  seconds. 

A  grid  search  using  GAUSS2  was  carried  out  for  selected 

values  of  a)  at  ^  =  0.5  and  T  =  1  second.  The  differences  for 
n 

e,  are  +0.4%  at  a)  =  0.1  and  -2.6%  at  w  = 

In  n 


96.0.  The 


. 
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FIGURE  5.14  PREDICTED  e±  FOR  PlCs  Of  G4(s)  WITH  T=1  SEC. 


NORMALIZED  NOISE  TRANSFER 
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NORMALIZED  NOISE  TRANSFERS  USING  PREDICTED 
PlCs  OF  G4(s)  WITH  T=1  SEC.  USING  GAUSS2 


PREDICTED  e1  FOR  PlCs  OF  G4(s)  WITH  T=5  SEC. 


FIGURE  5.16 
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normalized  noise  transfer  differences  are  -1.2%  and  -0.2% 
respectively.  This  agreement  is  considered  satisfactory. 


5-6  RESULTS  FOR  THE  FOUR  PLANTS  WHEN  A  RAMP 
OR  AN  ACCELERATION  INPUT  IS  USED 


Normally  when  a  controller  which  is  deadbeat  to  a 
unit  step  is  subjected  to  other  inputs,  errors  occur.  In 
particular,  the  response  to  a  ramp  or  an  acceleration  is 
very  poor.  The  unit  step  deadbeat  controller  exhibits  a 
constant  error  to  a  ramp  input  and  an  increasing  error  to  an 
acceleration  input. 

Comparisons  between  the  minimum  time  and  predicted 
PICs  for  the  four  plants  considered  are  presented  in  this 
section.  To  test  these  controllers,  the  mean  square  error 
between  the  input  and  output  was  measured  over  128  sampling 
periods.  The  error  for  the  predicted  PIC  was  then  normalized 
to  the  error  of  the  appropriate  minimum  time  controller. 

The  plants  tested  are 


1. 


2. 


3. 


4. 


(s)  with  T  = 
G2  (s)  with  T  = 
G^  (s )  with  T  = 
G^ (s)  with  T  = 


1  second  in  FIG. 
1  second  and  8  = 
1  second  and  6  = 
1  second  and  C  = 


5.17. 

1.0  in  FIG.  5.18. 
5.0  in  FIG.  5.19. 
0.5  in  FIG.  5.20. 


The  ramp  input  is 


ur(t) 


NT 


(5.24) 
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and  the  acceleration  input  is 

u  (t)  =  (NT)2  (5.25) 

a 

where  N  is  the  number  of  the  sampling  period  and 
0  ^  N  ^  128. 

The  graphs  shown  in  FIG.  5.17  to  5.19  are  plotted 
over  very  small  ranges  of  <*.  For  G^(s)  to  G^(s),  the 
improvement  of  the  PICs  over  the  minimum  time  controller 
occurs  only  in  this  region.  For  values  of  «  outside  these 
ranges,  the  error  of  the  PIC  was  larger.  In  particular  for 
«  >  1.0,  the  PIC  responds  more  slowly  than  the  minimum  time 
controller  and  therefore  the  normalized  error  rises  rapidly 
with  a. 

In  FIG.  5.20,  the  values  of  plotted  are  limited 
as  the  normalized  errors  are  greater  then  1.0  for  outside 
these  values.  This  is  similar  to  the  results  for  the  other 
plants  considered. 

For  the  four  plants  tested,  the  maximum  improvement 
in  the  error  is  significant.  The  maximum  improvement  is 
slightly  greater  for  the  acceleration  error.  The  range  of 
plant  parameters  for  which  this  improvement  is  obtained  are 
very  limited  in  each  case. 
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a 

FIGURE  5.17  NORMALIZED  RAMP  AND  ACCELERATION  ERRORS 

USING  PREDICTED  PlCs  OF  G-^s)  WITH  T=1  SEC. 

FIGURE  5.18  NORMALIZED  RAMP  AND  ACCELERATION  ERRORS 

USING  PREDICTED  PlCs  OF  G2(s)  WITH  6=1 
AND  T=1  SEC. 


NORMALIZED  PIC  ERROR 
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FIGURE  5.19  NORAMLIZED  RAMP  AND  ACCELERATION  ERRORS 

USING  PREDICTED  PlCs  OF  G-fs)  WITH  3=5 
AND  T=1  SEC. 


FIGURE  5.20  NORMALIZED  RAMP  AND  ACCELERATION  ERRORS 

USING  PREDICTED  PlCs  OF  G  (s)  WITH  C=0.5 
AND  T=1  SEC.  4 
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5-7  SOME  CONCLUSIONS  FOR  THE  PICs 

The  reduction  of  the  system  noise  transfer  varied 
from  approximately  25%  to  50%  when  the  predicted  PICs 
replaced  the  minimum  time  deadbeat  controllers  for  any  of 
the  plants  tested.  These  values  are  taken  from  the  normalized 
noise  transfer  curves,  obtained  using  the  noise  sources 
defined  in  Chapter  4,  for  the  various  plants  tested  using 
the  digital  simulation.  Similar  improvements  were  also 
observed  for  the  minimum  noise  transfer  controllers  of  the 
plants  tested  using  the  hybrid  simulation  in  Chapter  3. 

The  largest  improvements  occured  for  (s)  and  G^(s) 

which  have  4  term  PICs.  In  both  cases,  the  response  time  of 

the  PIC  was  1/3  longer  than  for  the  minimum  time  deadbeat. 

The  percentage  improvement  in  the  noise  transfer  was  generally 

numerically  greater  than  the  percentage  increase  in  the 

response  time  for  the  plants  tested.  For  some  sets  of 

parameters,  the  improvement  for  (s)  was  as  high  as  70%  and 

as  low  as  15%.  However  for  co  <  .5  and  co  >5,  the 

n  n 

improvements  fell  within  the  original  range  stated. 

For  G^(s)  and  G^(s),  the  percentage  increase  in  the 
response  time  caused  by  the  PICs  was  50%  as  n=3  for  the  PICs. 
In  these  cases,  the  reduction  of  the  noise  transfer  was  never 
numerically  >  50%. 

Therefore  increasing  the  response  time,  when  using 
the  predicted  PICs,  gave  the  best  noise  transfer  improvement 
for  higher  order  systems.  However,  the  improvement  for  all 
systems  tested  was  substantial  when  the  predicted  PICs  were 
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used.  Therefore  the  trade-off  of  increased  response  time 
and  greater  controller  complexity  for  improved  noise 
performance  is  justified  for  the  predicted  PICs. 

It  was  stated  earlier  that  the  input  noise  is 
assumed  to  be  white.  The  grid  searches  carried  out  for  each 
plant  used  the  non-white  noise  sources  defined  in  Chapter  4. 

From  these  tests,  two  conclusions  may  be  drawn. 

1.  For  all  plants  tested,  the  agreement  between  the 
predicted  and  "grid  searched"  PICs  was  generally 
very  close.  The  noise  transfer  differences 
were  always  <  3.7%  and  generally  were  <  2.0%. 

The  e^  differences  were  always  <  7.0%  for  any 
noise  source  used.  This  indicates  that  the 
predicted  PICs  are  at  least  a  good  approximation 
for  a  variety  of  non-white  noise  sources.  The 
results  from  the  hybrid  simulation  also  support 
this  statement. 

2.  The  agreement  between  the  predicted  and  "grid 

searched"  PICs  was  best  for  both  e^  and  the  noise 

transfers  when  the  GAUSS  generated  noise  sources 

were  used.  These  have  lower  deviations  from  a 

uniform  frequency  distribution  between  0  and 

go  /2  than  do  the  RANDU  generated  noise  sources, 
s 

This  supports  the  assumption  that  the  NCE 
criterion  is  derived  for  a  white  noise  input  as 
the  GAUSS  noise  sources  are  "whiter"  in  terms  of 
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their  frequency  distribution  than  are  the  RANDU 
noise  sources. 

For  the  plants  G-^(s),  G^fs)  and  G^(s),  define  the 
respective  dimensionless  constants  <*T,  “T  and  Coo^T  as  t. 

For  G>2  (s)  /  define  the  two  dimensionless  constants  <*T  and 
$T  as  x^  and  These  constants  are  analogous  in  form  to 

time  constants. 

In  this  chapter,  all  noise  transfer  curves  were 
plotted  for  a  sampling  rate  of  1  second.  All  these  curves 
had  one  common  tendency.  When  t  was  very  large  or  very 
small,  the  rate  of  change  of  the  noise  transfer  curve  for 
changes  in  t  was  very  slow.  Between  these  two  regions 
there  was  a  transition  zone.  The  transition  zone  extended 

(\j  f\f 

from  x  ^  0 .  IT  to  x  'v  10T  in  all  cases. 

One  reason  for  this  transition  zone  is  the  interaction 
of  the  values  in  t.  For  x  <  0.1T,  the  dominant  terms  are  « , 
or  3  depending  on  the  plant  under  consideration.  For 
x  >  10T,  T  is  the  dominant  term.  In  the  transition  zone, 
there  is  no  clearly  dominant  term. 

Two  constants  x^  and  x^  exist  for  G ^  (s) .  Either  one 
may  cause  a  transition  zone  to  exist  if  the  other  is  held 
constant . 

Within  this  transition  zone,  the  absolute  maximum  and 
the  absolute  minimum  noise  reductions  occurred  for  all  plants 
tested.  For  G^(s),  the  maximum  reduction  occurred  at  x  =  0.7 
and  the  minimum  reduction  occurred  at  x=3 . 5 .  The  appropriate 
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values  of  t  may  also  be  found  for  the  other  plants  tested. 

Therefore  a  set  of  plant  parameters  exists  which  give 
the  best  possible  noise  improvement  for  a  specific  plant 
type.  If  the  plant  parameters  or  the  system  sampling  rate 
may  be  altered  to  any  extent,  the  constant  t  can  be  adjusted 
to  give  the  best  noise  reduction. 

The  sensitivity  of  a  deadbeat  controller  is  effectively 
infinity  as  the  response  to  a  unit  step  is  either  deadbeat  or 
it  is  not.  However  if  a  small  parameter  variation  occurred 
in  the  plant,  the  difference  between  the  ideal  deadbeat 
response  and  the  perturbed  response  was  generally  not  very 
large.  Both  the  minimum  time  deadbeat  controllers  and  the 
predicted  PICs  were  tested  for  this  error. 

For  the  four  plants,  one  of  the  parameters  was  varied 
by  ±  5%  and  the  resultant  errors  in  the  response  to  a  unit 
step  were  evaluated  over  128  sampling  periods.  For  G^(s), 

G2(s)  and  G^ (s) ,  this  parameter  was  a.  For  G^(s),  was 
varied.  In  all  cases,  the  differences  between  the  errors 
measured  for  the  minimum  time  controllers  and  the  predicted 
PICs  were  less  than  2%.  Depending  on  the  parameters,  the 
error  could  be  lower  for  either  type  of  controller.  Therefore 
the  predicted  PICs  caused  no  degradation  in  this  respect 
when  compared  to  the  corresponding  minimum  time  deadbeat 


controller . 
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CHAPTER  VI 


CONTROLLERS  EXTENDED  TWO  TERMS 

In  this  chapter,  the  results  of  the  purely  digital 
simulation  of  the  controllers  extended  by  two  terms  are 
presented.  To  simplify  the  presentation,  a  controller 
extended  by  two  terms  is  referred  to  as  a  P2C.  This 
abreviation  stands  for  Plus  Two  Terms  Controller. 

The  algorithim  used  to  calculate  the  P2Cs  is  derived 
in  the  first  section.  This  derivation  is  based  on  the  design 
criterion  NCE.  In  the  following  sections,  the  results 
obtained  from  the  digital  simulation  of  the  four  plants  listed 
in  Chapter  1  are  presented.  These  results  are  presented  in 
the  same  format  used  in  Chapter  5.  To  compare  the  PICs  and 
the  P2Cs,  a  graph  showing  the  normalized  noise  transfers 
obtained  using  both  types  of  predicted  controllers  is  given 
for  each  plant  tested.  The  final  section  of  the  chapter 
presents  some  conclusions  derived  from  the  test  results. 


6-1  ALGORITHIM  USED  IN  THE  DIGITAL  SIMULATION 
TO  CALCULATE  THE  PREDICTED  P2Cs 

In  this  section,  the  algorithim  used  to  calculate 
the  predicted  P2Cs  is  derived  from  NCE.  It  was  shown  in 
Section  5-2,  that  one  extra  term  gave  one  additional  degree 
of  freedom  in  the  calculation  of  the  D(z).  This  degree  of 
freedom  was  represented  by  the  choice  of  e^. 
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Similarly  adding  two  extra  terms  to  the  controller 
corresponds  to  two  extra  degrees  of  freedom  in  the 
calculation.  The  choices  of  e^  and  e 2  represent  these  two 
degrees  of  freedom. 

The  controllers  are  calculated  by  the  method  derived 
in  Section  2-4.  Expressing  the  remaining  coefficients  of 
the  P2C  as  functions  of  e^  and  e 2  yields: 


h 

h 

h 

h 

h 

e 

e 


0 

1 

2 

3 

4 

3 

4 


alel  +  a2e2  + 


a4el  +  a5e2  + 


a7el  +  a8e2  + 
a10el  +  alle2 
a13el  +  a14e2 


a6 

a9 

+  a12 
+  a15 


blel  +  b2S2  +  b3 
b4el  +  b5S2  +  b6 


(6.1) 

(6.2) 

(6.3) 

(6.4) 

(6.5) 

(6.6) 

(6.7) 


where  the  order  of  the  P2C  is  n  ^  4.  These 
coefficients  are  given  in  Table  6.1  for  the  four  plants 
tested.  The  coefficients  are  defined  in  terms  of  the 
constants  derived  in  Chapter  1  for  each  plant. 

For  controllers  with  n  <  4,  the  form  of  the  criterion 


NCE  is 
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NCE  h0  +  ^hl  elh0^  +  ^h2  elhl  e2h0  + 


el  ho^  +  (h3  elh2  e2hl  e3h0  + 


2  3  2 

el  hl  “  el  h0  +  2ele2h0^  +  ^h4  “  elh3 


e2h2  e3hl  e4h0  +  el  h2  +  2ele2hl  + 


2ele3h0  +  e2  h0  3el  e2h0  el  hl  +el  V 


(6.8) 


Since  the  coefficients  defined  in  (6.1)  to  (6.7)  are 
functions  of  e^  and  e^r  NCE  may  be  written  in  the  same  form. 
Define  NCE£  as  the  criterion  for  the  system  with  two  extra 
sampling  periods.  Rewriting  NCE  in  this  form  now  yields: 


n 

NCE9  =  f  (e  ,  e„)  =  I  NCE 

z  1  i=0 


where 


NCE 


20 


(6.9) 


(6.10) 


NCE24  (h4  elh3  e2h2  e3hl  ~  e4h0  + 


el  h2  +  2ele2hl  +  2ele3h0  +  e2  h0 


,  2  .  3,  4,  >2 

3el  e2h0  el  hl  +  el  h0^ 


(6.11) 


eisdw 
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The  substitution  of  (6.1)  to  (6.7)  into  (6.8) 
yields  NCE^  which  contains  products  of  the  powers  of  e^  and 
e 2*  Thus  it  is  not  possible  to  find  the  values  of  e^  and 
which  correspond  to  the  minimum  of  NCE2  by  simple 
differentiation  as  in  Chapter  5. 

An  alternate  method  using  the  partial  differentiation 
of  NCE2  with  respect  to  e^  for  fixed  values  of  e^  is  derived 
in  the  remainder  of  this  section. 

Substituting  (6.1)  to  (6.7)  into  (6.8)  yields  in 
simplified  form: 


NCE20  =  a22e22  +  2a2b01e2  +  b012 


NCE21  bll  e2  +  2bllb1262  +  b12 


NCE22  a2  e2  2a2b21e2  +  ^b21 


2a2b22)e2  +  2b21b22S2  +  b22 


NCE23  “  b31  e2  +  2b31b32e2  +  (b32  + 


2b31b33)e2  +  2b32b33e2  +  b33 


NCE2  4  "  b41  e2  +  2b41b42e2  +  (b42  + 


2b41b43)e2  +  2b42b43e2  +  b43 


(6.12) 

(6.13) 

(6.14) 

(6.15) 

(6.16) 
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where 


b01 

alel  +  a3 

bll  1 

a5  "  a2el 

b12 

=  -aiei2  +  (a4  -  a3)ei  +  ag 

b2 1 

2 

=  (ag  -  a9)  -  (a1  +  a5)e1  +  a2e1 

b22 

2  3 

~  a  ^  “I-  ( a  —  a^ )  e  ^  (a3  —  a4 )  ei  a  ^  e 

b31 

2a-^e^  -  a2b2  ~  a^ 

b32 

3  2 

=  “a2ei  +  ^ai  +  a5^el  +  ^-a4  “  ^alb2  + 

b33 

2a3  '  a8)el  +  (all  "  a6  "  2a2b3) 

4  3  2 

alal  "**  (a4  ~  a.3 ) ei  (a^  ~  3^  ~  a^b^)©^ 

(a10  a9  2alb3^el  +  ^ a12  “  a3b3^ 

b41  = 

2 

=  ^a^^^  +  (a1  +  2a2b2  +  2a5)e1  +  (-ag  -  a2b5  - 

b42 

a5b2  +  a3) 

=  a2e].4  ~  ^a5  +  3a1)e13  +  (4a3b2  -  3a3  + 

2 

2a4  +  ag )  ei  +  ("an  "  a7  ”  2a4b2  ~  2alb5  + 

+  2a6  +  4a2b3)e;L  +  (a14  -  ag  -  2a5b3  -  2a2bg) 

. 
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b4  3 

:  a1e15  +  (a3 

i 

pj 

CD 

I-1 

+  (2a^b1  +  a^  -  a6 

)e13  + 

^4alb3  +  a9 

-  axb4  - 

a10  -  a4bl)el2  + 

(a13  "  a12 

-  2a4b3  - 

2alb6  +  ^a3b3^el  + 

(a15  "  a6b3 

a3b6* 

(6.17) 

Collecting  (6. 

12)  to  (6.16) 

yields  the  simplified 

expression : 

NCE0 

n 

=  I  NCE„ . 

i-o  21 

=  D4e24 

3  2 

+  D_e0  +  D0e~  + 

j  Z  A  A 

°le2  +  D0 

(6.18) 

where  2  <  n  ^ 

4  and 

D. 

l 

fi(e1)  ;  i 

=  0,  1,  . 

4 

(6.19) 

„  _  <  <  , 

For  0  n  n  ^  1 

,  NCE  cannot 

be  calculated  as  no  P2C  can 

exist  for  these  values. 

When  Dq  to  are  evaluated,  then  (6.18)  is 
expressed  only  as  a  function  of  e^.  If  the  values  of  e^ 
are  known,  then  Dq  to  may  be  calculated.  One  way  of 
obtaining  these  values  of  e^  is  by  iteration  over  0  ^  e^ 
\  1.0  which  is  required  for  deadbeat. 

For  a  particular  value  of  e^,  equation  (6.18)  may 
be  rewritten  in  simplified  form  as: 


* 
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NCE2 


[f  (e2) ] 


(6.20) 


Setting  the  partial  derivative  of  (6.20)  with  respect  to  e2 
to  zero  yields: 


[nce2] 


[f  (e2) ] 


Substituting  (6.18)  into  (6.21)  yields: 


(6.21) 


4(D4 

)e23  +  3(D3 

)e22  +  2(D2 

ei 

e  -t 

>  e2  + 


(D. 


)  =  0 


(6.22) 


From  (6.22),  solving  for  the  value  of  e2  which  satisfies  the 
deadbeat  requirements: 


0  ^  e2  ^  1.0  (6.23a) 

e2  '  e1  (6.23b) 

yields  the  value  of  e2  which  corresponds  to  the  minimum 
value  of  NCE2  for  a  particular  value  of  e^.  The  criterion 
NCE2  may  be  evaluated  by  substituting  the  value  of  e^  and 
the  corresponding  value  of  e2  into  (6.18)  and  (6.19). 

From  the  values  of  NCE2  calculated  for  all  e^  values 
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obtained  by  iteration,  the  global  minimum  value  of  NCE2  may 
be  selected.  The  values  of  e^  and  e^  which  correspond  to 
this  value  of  NCE^  are  then  the  predicted  values  for  the 
plant  under  consideration. 

The  method  of  calculation  may  be  outlined  as  follows. 

1.  Write  NCE2  entirely  as  a  polynominal  in  e2 
where  the  coefficients  are  functions  of  e^. 

2.  Iterate  e^  over  the  range  0  to  1.0  to  calculate 
these  coefficients. 

3.  For  each  e^  value,  take  the  partial  derivative 
of  NCE^  with  respect  to  e^  and  set  it  to  zero. 
Solve  for  the  real  root  of  e2  which  satisfies 
(6.23) . 

4.  Evaluate  NCE2  at  each  value  of  e^  and  the 
corresponding  value  of  e2*  The  value  of  NCE2 
which  is  the  global  minimum  yields  the  e^  and 
e2  values  of  the  predicted  P2C. 

6-2  RESULTS  FOR  P2Cs  OF  G^s)  =  1/s  (s  +  <*) 

This  plant  was  tested  with  its  P2Cs  using  the  same 
test  procedure  as  used  for  the  PICs  in  Chapter  5. 

The  sampling  rate  was  set  at  T  =  1  second  and  a 
was  varied  from  0.125  to  32.0.  The  predicted  values  of  e^ 
and  e2  are  shown  in  FIG.  6.1  plotted  versus  a. 

FIG.  6.2  is  a  plot  of  the  normalized  noise  transfers 
of  the  systems  containing  the  predicted  P2Cs  versus  *.  To 
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a 

FIGURE  6.1  PREDICTED  e,  AND  e„  FOR  P2Cs  OF  G, (s)  WITH 

T=1  SEC.  1 

FIGURE  6.2  NORMALIZED  NOISE  TRANSFERS  USING  PREDICTED 

PICs  AND  P2Cs  OF  G1(s)  WITH  T=1  SEC.  USING 

GAUSS 2 
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FIGURE  6.3  PREDICTED  e1  AND  e2  AND  GRID  SEARCHED  e1  AND 

e2,  USING  GAUSS 2 ,  FOR  P2Cs  OF  G-^s)  WITH  T=1  SEC. 

FIGURE  6.4  NORMALIZED  NOISE  TRANSFERS  USING  PREDICTED 

AND  GRID  SEARCHED  P2Cs  OF  G, (s)  WITH  T=1  SEC. 
USING  GAUSS2  1 
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obtain  the  comparison  between  the  PlCs  and  P2Cs,  the 

normalized  noise  transfers  obtained  using  the  predicted 

PlCs  with  T  =  1  second  are  also  plotted  on  FIG.  6.2. 

Next  a  grid  search  was  performed  using  the  noise 

source  GAUSS2.  The  comparison  between  the  predicted  values 

of  e^  and  e^  and  the  grid  searched  values  is  shown  in 

FIG.  6.3.  FIG.  6.4  shows  the  comparison  between  the 

normalized  noise  transfers  for  the  systems  containing  the 

grid  searched  and  predicted  P2Cs. 

The  improvement  of  the  P2Cs  over  the  PlCs  varies 

between  7%  and  12%  as  calculated  by  GAUSS2  when  referred  to 

the  minimum  time  deadbeat.  From  FIG.  6.2,  the  greatest 
.  v 

improvements  occur  at  «  ^  0.5  and  for  «  >  10.0. 

The  predicted  and  grid  searched  results  in  this 
chapter  are  compared  in  the  same  way  as  in  Chapter  5.  The 
e^  differences  are  -7.0%  at  a  =  0.125  and  -5.3%  at  «  =  32.0. 
Respectively,  the  e 2  differences  are  -27%  and  +17%.  The 
corresponding  transfer  differences  are  -20%  and  -5.3%. 

These  noise  transfer  results  exhibit  the  same 
tendencies  as  the  transfers  for  the  PlCs.  However,  the 
variations  involved  are  much  larger. 

6-3  RESULTS  FOR  P2Cs  OF  G2(s)  =  1/s  (s+«)  (s+3) 

For  this  plant,  the  values  tested  were  3  =  1.0  and 
T  =  1  second  with  varied  from  0.2  to  32.0. 

FIG.  6.5  gives  the  predicted  values  of  e^  and  e2  for 
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this  plant.  FIG.  6.6  shows  the  comparison  between  the 
normalized  noise  transfers  obtained  using  the  predicted 
PlCs  and  P2Cs  for  a  GAUSS2  input.  In  FIG.  6.7,  the  grid 
search  calculated  values  of  e^  and  e^  using  GAUSS2  are 
compared  to  the  predicted  values.  The  comparison  of  the 
normalized  noise  transfers  for  the  systems  using  the 
predicted  and  grid  searched  P2Cs  is  shown  in  FIG.  6.8. 

A  difficulty  in  calculating  the  predicted  values  of 
e^  and  e 2  for  “  ^  1.0  was  encountered.  The  solution,  as 
previously  outlined,  gave  erroneous  results  for  e2»  To 
overcome  this  problem,  the  values  of  e^  and  e^  were  both 
iterated  to  find  the  minimum  value  of  NCE2.  This  value  then 
yielded  the  predicted  e^  and  e 2  values  plotted.  For  values 
of  «  >  1.0,  no  difficulty  was  encountered.  To  check  the 
results  for  a  >  1.0,  the  values  of  e^  and  e2  were  iterated 
and  the  global  minimum  value  of  NCE2  calculated.  For  the 
values  of  a  >  1.0  checked,  the  results  of  both  methods 
agreed. 

The  improvement  in  the  normalized  noise  transfers 
when  using  the  P2Cs  in  place  of  the  PlCs  is  from  7%  to  13% 
as  calculated  using  GAUSS2.  This  is  taken  from  FIG.  6.6. 

The  e^  differences  from  FIG.  6.7  are  -2.2%  at  a  =  0.4 
and  -8.3%  at  a  =  32.0.  For  e2,  these  differences  are  -7% 
and  -15.2%  respectively.  The  noise  transfer  differences  are 
-6.8%  and  -3.6%  respectively. 

For  this  plant,  the  differences  between  the  noise 
transfers  are  smaller  than  the  error  term  differences.  This 
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a 

FIGURE  6.5  PREDICTED  e  AND  e2  FOR  P2Cs  OF  G2(s)  WITH 

3=1  AND  T=1  SEC. 

FIGURE  6.6  NORMALIZED  NOISE  TRANSFERS  USING  PREDICTED 

PICs  AND  P2Cs  OF  G?(s)  WITH  3=1  AND  T=1  SEC. 
USING  GAUSS2 
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PREDICTED 


_ i _ i _ i _ 

0  1  10  10  0  1000 

a 

6.7  PREDICTED  e_L  AND  e2  AND  GRID  SEARCHED  e1  AND 
e2,  USING  GAUSS 2 ,  FOR  P2Cs  OF  G2(s)  WITH  6=1 
AND  T=1  SEC. 

6.8  NORMALIZED  NOISE  TRANSFERS  USING  PREDICTED 
AND  GRID  SEARCHED  P2Cs  OF  G2(s)  WITH  6=1  AND 
T=1  SEC.  USING  GAUSS 2 
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indicates  that  the  rate  of  change  of  the  noise  transfer 
with  e^  and  e 2  is  rather  slow. 


6-4  RESULTS  FOR  P2Cs  of  G3(s)  =  (s+3)/s(s+«) 


The  values  considered  in  this  case  were  3  =  5.0  and 
T  =  1  second  with  varied  from  0.25  to  64.0.  FIG.  6.9 
illustrates  the  predicted  values  of  e^  and  e^.  The  normalized 
noise  transfers  obtained  using  the  predicted  PlCs  and  P2Cs  as 
calculated  for  GAUSS2  are  shown  in  FIG.  6.10.  FIG.  6.11  and 
FIG.  6.12  compare  the  predicted  and  grid  searched  P2Cs  error 
terms  and  the  normalized  noise  transfers  as  calculated  for 
GAUSS2  respectively. 

The  actual  improvement  of  the  P2Cs  over  the  PIC  varies 
from  7.0%  at  «  =  64.0  to  11.5%  at  «  =  0.5.  These  figures  are 
taken  from  FIG.  6.10. 

From  FIG.  6.11,  the  e^  differences  are  -13.4%  at 
«  =  0.25  and  -5.5%  at  «  =  32.0.  Similarly  the  e 2  differences 
are  -30%  and  +  28.1%  respectively.  The  corresponding 
differences  for  the  noise  transfers  are  -10.0%  and  -3.2%. 

Again,  there  are  large  differences  for  the  error 
terms.  However  the  differences  for  the  normalized  noise 
transfers  are  substantially  smaller. 

2  2 

6-5  RESULTS  FOR  P2Cs  OF  G4(s)  =  l/s(s  +2CwnS  +  ) 

This  plant  was  tested  for  £  =  0.5  and  T  =  1  second 

with  a)  varied  from  0.25  to  64.0. 
n 
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a 

FIGURE  6.9  PREDICTED  e±  AND  e2  FOR  P2Cs  of  G3(s)  WITH 

3=5  AND  T=1  SEC. 

FIGURE  6.10  NORMALIZED  NOISE  TRANSFERS  USING  PREDICTED 

PlCs  AND  P2Cs  OF  G3(s)  WITH  3=5  AND  T=1  SEC. 


NORMALIZED  NOISE  TRANSFER  D(z)  ERROR  TERMS 


135 


FIGURE  6.11  PREDICTED  e1  AND  e2  AND  GRID  SEARCHED  e1  AND 

e2,  USING  GAUSS 2 ,  FOR  P2Cs  OF  G3(s)  WITH  3=5 
AND  T=1  SEC. 

FIGURE  6.12  NORMALIZED  NOISE  TRANSFERS  USING  PREDICTED  AND 

GRID  SEARCHED  P2Cs  OF  G  (s)  WITH  6=5  AND  T=1 
SEC.  USING  GAUSS2  J 
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6.13  PREDICTED  e1  AND  e2  FOR  P2Cs  OF  G4 (s)  WITH 
£=0.5  AND  T=1  SEC. 

6.14  NORMALIZED  NOISE  TRANSFERS  USING  PREDICTED 
PlCs  AND  P2Cs  OF  G.(s)  WITH  £=0.5  AND  T=1 
SEC.  USING  GAUSS 2 4 
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FIGURE  6.15  PREDICTED  e ±  AND  e2  AND  GRID  SEARCHED  e1  AND 

e2,  USING  GAUSS 2 ,  FOR  P2Cs  OF  G4(s)  WITH 
c=0 . 5  AND  T=1  SEC. 

FIGURE  6.16  NORMALIZED  NOISE  TRANSFERS  USING  PREDICTED 

AND  GRID  SEARCHED  P2Cs  OF  G4  (s)  WITH  C=0.5 

AND  T=1  SEC.  USING  GAUSS2 
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As  for  G2  (s)  ,  the  solution  was  not  obtainable  by  use 

of  the  calculation  algorithim  for  some  values.  When  u)  <2, 

n 

the  predicted  values  of  e^  and  e2  were  calculated  by  iterating 
until  the  minimum  value  of  NCE2  was  found.  For  values  of 
a)  >  2,  no  such  difficulty  existed. 

FIG.  6.13  shows  the  predicted  values  of  e^  and  e^. 

The  normalized  noise  transfers  obtained  using  the  predicted 
PICs  and  P2Cs  when  subjected  to  GAUSS2  are  compared  in  FIG. 
6.14.  A  comparison  of  the  D(z)  error  terms  for  the  predicted 
and  grid  searched  P2Cs  is  given  in  FIG.  6.15.  The  final 
graph,  FIG.  6.16  compares  the  normalized  noise  transfers 
obtained  using  the  predicted  and  grid  searched  P2Cs  using 
GAUSS 2 . 

The  improvement  in  the  normalized  noise  transfers 
when  using  the  predicted  P2Cs  instead  of  the  predicted  PICs 

f\j 

varies  from  essentially  zero  at  con  ^  1.5  to  approximately 
13%  at  03^  =  0.25.  For  >  10.0,  the  improvement  is  greater 
than  10%. 

The  differences  in  e,  are  -1%  at  (1)  =0.25  and  -12% 

1  n 

at  u)  =64.0.  The  differences  for  e~  are  -5.9%  and  +  11% 
n  2 

respectively  as  calculated  from  FIG.  6.15.  From  FIG.  6.16, 

the  maximum  noise  transfer  difference  is  -14.7%  at  w  =0.5 

n 

v 

and  the  minimum  difference  is  zero  for  w  v  1.25.  For 

n 

00  >4,  the  difference  is  no  greater  then  -3%. 

n 
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6-6  SOME  CONCLUSIONS  FOR  THE  P2Cs 

The  reduction  of  the  system  noise  transfer  varied 
from  25%  to  70%  when  the  predicted  P2Cs  replaced  the  minimum 
time  deadbeat  controllers  for  any  of  the  plants  tested.  These 
percentages  were  calculated  using  GAUSS2  as  the  noise  source. 

The  largest  improvements  occurred  for  (s)  and  G^  (s) 
which  had  five  term  P2Cs.  This  represented  a  2/3  increase  in  the 
response  time  when  compared  to  the  minimum  time  deadbeat 
controller.  The  ranges  of  improvements  were  40%  to  66%  for 
G2(s)  and  25%  to  70%  for  G^(s). 

When  the  larger  noise  transfer  reductions  are 
considered,  the  trade-off  of  response  time  would  generally 
be  justified.  However  for  some  regions,  the  numerical 
value  of  the  percentage  noise  reduction  achieved  was  much  less 
than  the  numerical  value  of  the  percentage  increase  in  the 
response  time.  In  these  regions,  the  decision  must  be  based 
on  the  importance  of  the  absolute  noise  reduction. 

For  G^(s)  and  G^(s) ,  the  increase  in  the  response 
time  was  100%.  The  ranges  of  the  noise  reductions  were  35% 
to  45%  for  G^(s)  and  40%  to  55%  for  G^[s)  .  In  both  cases, 
the  numerical  value  of  the  percentage  noise  reduction  was 
much  less  than  the  numerical  value  of  the  percentage  increase 
in  the  response  time.  Therefore  the  decision  to  use  these 
P2Cs  must  be  based  on  the  importance  of  the  absolute  noise 
reduction . 

Again  for  the  predicted  P2Cs,  a  transition  region 
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in  the  noise  transfer  curves  existed  for  the  same  values 
of  t  as  found  for  the  PICs.  The  shape  of  the  noise  transfer 
curves  using  the  predicted  P2Cs  was  virtually  the  same  as  for 
the  corresponding  predicted  PICs  for  all  plants  tested.  The 
absolute  values  of  the  noise  reduction  were  the  only 
differences.  The  reasons  for  the  shape  of  these  curves  are 
the  same  as  those  given  in  Section  5-7. 

The  agreement  between  the  predicted  and  "grid 
searched"  P2Cs  was  not  as  good  as  the  agreement  obtained  for 
the  PICs.  Generally  the  agreement  between  the  e^  values  was 
better  than  the  agreement  between  the  e 2  values. 

The  maximum  difference  for  the  noise  transfers  was 
20%  but  the  difference  generally  was  <  10%.  For  t  <  0.5T, 
the  differences  were  the  largest.  These  differences  dropped 
to  <  5%  for  x  >  T  for  all  plants  tested. 

The  differences  between  the  e^  values  were  generally 
comparable  to  the  noise  transfer  differences.  However,  the 
differences  between  the  e 2  values  were  generally  numerically 
larger  than  the  noise  transfer  differences. 

Since  GAUSS2  is  not  a  pure  white  noise  source,  it 
has  some  correlation  with  the  step  input.  The  P2C  acts  on 
the  system  for  a  longer  period  of  time  than  does  the 
corresponding  PIC.  Therefore  any  differences  due  to  the 
non-white  noise  source  would  be  magnified  by  the  P2C .  This 
is  the  reason  that  the  differences  between  the  predicted 
and  "grid  searched"  P2Cs  are  larger  than  for  the  PICs. 
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For  large  values  of  t,  the  system  response  occurs 
mostly  in  the  first  part  of  the  controller  cycle.  Therefore 
changes  in  the  later  terms  of  the  controller  have  a  lesser 
effect  than  when  t  is  small.  This  accounts  for  the  better 
agreement  between  the  "grid  searched"  and  predicted  P2Cs  for 
larger  values  of  t.  For  t  >  T,  the  predicted  P2Cs  appear  to 
be  an  acceptable  approximation  of  the  minimum  noise  transfer 
controllers  for  the  non-white  noise  sources  used.  For  t  <  T, 
only  G2(s)  exhibits  any  continuing  agreement  between  the 
predicted  and  "grid  searched"  P2Cs.  The  agreement  for  the 
P2Cs  is  never  as  good  as  the  agreement  for  the  PICs  for  any 
of  the  plants  tested. 

Again  there  were  absolute  maximum  and  absolute 
minimum  noise  reductions  in  the  noise  transfer  curve  for  all 
plants  tested.  The  same  discussion  used  in  Section  5-7  may 
be  applied  in  this  case. 


CHAPTER  VII 


OTHER  EXTENDED  CONTROLLERS 

In  the  first  section,  the  predicted  controllers 
extended  by  3  terms  for  G^(s)  are  presented.  No  grid 
search  calculations  are  presented  for  this  case.  This 
controller  is  abreviated  as  P3C  which  stands  for  Plus  Three 
Extra  Terms  Controller. 

In  the  second  section,  some  predicted  controllers 
using  a  reduced  sampling  period  are  studied.  Some  results 
for  the  straight  line  approximations  to  the  predicted  PlCs 
are  presented  in  the  third  section. 

In  each  section,  some  conclusions  are  presented  for 
the  particular  situation  under  consideration. 

7-1  PREDICTED  P3Cs  OF  G^s) 

For  these  calculations,  T  =  1  second  was  the  sampling 
rate  and  was  varied  from  0.25  to  6  4.0. 

Adding  three  extra  terms  to  the  controller  results  in 
three  degrees  of  freedom  in  its  calculation.  This  conclusion 
is  reached  in  the  same  way  as  for  the  PlCs  and  P2Cs.  These 
three  degrees  of  freedom  may  be  represented  by  e^,  e 2  and  e^ 
in  the  denominator  of  the  controller. 

To  calculate  the  predicted  P3Cs,  the  expression  for 
NCE  when  n  =  5  was  written  from  (3.19).  Using  a  method 
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analogous  to  that  used  in  Chapters  5  and  6,  the  coefficients 
11q,  h^,  h 2/  h^ ,  h^  and  e^  were  written  as  functions  of  e^, 
e^  and  e^.  Using  these  expressions,  (3.19)  was  rewritten 
as  a  function  of  e^,  e 2  and  e^. 

By  iterating  e^,  e 2  and  e^  in  that  order  a  grid 
search  for  the  global  minimum  value  of  NCE  at  each  value  of 
a  was  performed.  This  approach  was  used  rather  than  the 
calculation  algorithim  used  in  Chapters  5  and  6  as  only  a 
few  results  were  taken.  The  programming  complexity 
associated  with  the  algorithim  was  thus  replaced  by  a 
simpler  iteration  process.  The  values  obtained  for  e^,  e ^ 
and  e^  are  given  only  to  2  significant  figures.  Once  the 
predicted  values  were  known,  the  resultant  controllers  were 
tested  using  the  GAUSS2  noise  input. 

The  predicted  values  of  e^,  e 2  and  e^  are  plotted  in 
FIG.  7.1.  A  comparison  of  the  normalized  noise  transfers  of 
the  systems  containing  the  predicted  PICs,  P2Cs  and  P3Cs  is 
shown  in  FIG.  7.2.  GAUSS 2  was  the  noise  source  used  in  all 
cases . 

Only  the  predicted  P3Cs  were  calculated  for  this 
plant  as  these  would  indicate  the  improvement  obtained  when 
the  P2Cs  are  replaced  by  the  P3Cs.  Also  a  grid  search  using 
GAUSS2  is  very  costly  in  terms  of  computing  time  and  therefore 
was  not  performed. 

The  improvements  in  the  normalized  noise  transfer 
when  using  the  P2Cs  rather  than  the  PICs  varied  from  6%  to 
13%  as  taken  from  FIG.  7.2.  The  corresponding  improvements 
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when  the  P3Cs  replace  the  P2Cs  are  0.5%  to  6%.  These 
percentages  are  measured  with  respect  to  the  corresponding 
minimum  time  controller. 

For  most  situations,  the  advantages  of  the  lowered 
noise  transfers  obtained  when  using  the  P3C  are  outweighed 
by  the  increased  complexity,  slow  response  time  to  a  stop 
input  and  the  cost  of  realization  of  the  P3Cs.  Therefore 
in  general  the  P3Cs  is  not  useful. 

However  for  some  particular  situations,  where  long 
term  noise  performance  is  the  most  important  requirement, 
the  use  of  these  P3Cs  may  be  warranted.  From  FIG.  7.2,  the 
noise  improvement  is  >  3%  for  >  1.0.  This  is  the  region 
where  the  largest  noise  improvement  occurs  as  compared  to 
the  P2C . 

7-2  PlCs  WITH  A  REDUCED  SAMPLING  PERIOD 

The  extended  controllers  considered  so  far  have  used 
the  same  sampling  rate  as  the  minimum  time  deadbeat  controller. 
The  case  will  now  be  considered  where  the  sampling  rate  used 
with  the  PIC  is  such  that  the  total  time  to  deadbeat  for  both 
the  minimum  time  controller  and  the  PIC  is  the  same. 

The  total  response  time  of  the  minimum  controller  is 
nT.  Therefore  the  sampling  period  of  the  PIC  must  be  reduced 
to  (^py)T,  where  n  is  the  order  of  the  minimum  controller. 

In  this  way,  the  total  response  time  to  deadbeat  is  constant 
for  both  controller  types. 
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FIGURE  7.1  PREDICTED  e2  and  e3  FOR  P3Cs  OF  G1(s) 

WITH  T=1  SEC. 

FIGURE  7.2  NORMALIZED  NOISE  TRANSFERS  USING  PREDICTED 

PlCs ,  P 2Cs  AND  P 3Cs  OF  G1(s)  WITH  T=1  SEC. 
USING  GAUSS2 
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In  some  situations,  it  may  be  desirable  that  the 
total  response  time  of  the  system  not  be  altered.  However 
some  noise  reduction  may  be  needed.  By  reducing  the  sampling 
period,  an  extended  controller  can  be  substituted  for  a 
minimum  time  controller  to  improve  the  noise  performance.  In 
this  section,  results  of  some  investigations  into  this  method 
are  presented. 

The  plants  considered  are  G-^(s)  and  G^  (s)  where  the 
sampling  period  with  the  minimum  time  controller  is  1  second. 
For  G^(s)  and  G£ (s)  respectively,  the  sampling  periods  with 
the  PICs  are  2/3  and  3/4  second. 

The  noise  transfers  obtained  using  the  predicted  PlCs 
with  T  =  1  second  and  the  appropriate  increased  sampling  rate 
are  calculated.  To  allow  direct  comparison,  both  types  of 
PICs  are  normalized  to  the  noise  transfer  obtained  using  the 
minimum  deadbeat  controller  with  T  =  1  second.  By  performing 
the  comparison  in  this  fashion,  the  effect  of  increasing  the 
sampling  rate  should  be  evident. 

To  allow  these  plants  to  be  tested  using  the  digital 
simulation,  appropriate  noise  sources  must  be  defined.  In 
practice,  the  minimum  time  and  extended  controllers  do  not 
"see"  exactly  the  same  noise  values  at  the  respective 
sampling  instants.  This  is  due  to  the  change  in  the  sampling 
rate.  For  G^(s),  the  sampling  instants  for  the  minimum 
time  and  extended  controllers  will  only  coincide  every  2T 
periods  and  similarly  every  3T  periods  for  G^  (s) . 
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A  sequence  of  random  numbers  is  generated  from  which 
the  appropriate  noises  sources  are  taken.  The  total  sequence 
generated  is  GAUSS 2A. 

For  G^(s),  a  random  number  was  generated  every  1/3 
second  in  the  simulation.  Therefore  every  third  number  was 
chosen  to  simulate  the  noise  source  when  T  =  1  second. 
Similarly  every  second  number  was  chosen  to  simulate  the 
noise  source  when  T  =  2/3  second.  The  experimental  run  in 
both  cases  was  256  seconds.  The  results  for  the  predicted 
PICs  using  both  sampling  rates  are  shown  in  FIG.  7.3. 

For  G2(s),  a  random  number  was  generated  every  1/4 
second  in  the  simulation.  When  T  =  1  second,  every  fourth 
number  was  chosen.  When  T  =  3/4  second,  every  third  number 
was  chosen.  Again  the  length  of  the  experimental  run  was 
256  seconds.  These  results  are  shown  in  FIG.  7.4. 

The  noise  sources  used  in  these  simulations  were 
called  GAUSS 2AM . 

From  FIG.  7.3  and  FIG.  7.4,  the  noise  reduction 
obtained  when  using  the  increased  sampling  rate  predicted 
PICs  is  not  as  great  as  that  obtained  when  using  the  normal 
predicted  PICs.  For  both  plants,  this  difference  is  largest 
for  small  values  of  a.  The  results  for  the  two  types  of 
PICs  agree  more  closely  as  a  increases.  For  «  >  3  the 
increased  sampling  rate  predicted  PIC  is  slightly  better  for 
G^  (s ) . 

A  given  plant  requires  that  the  controllers  have 
higher  equivalent  gains  during  the  sampling  period  when  the 
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FIGURE  7.3  NORMALIZED  NOISE  TRANSFERS  USING  PREDICTED 

PlCs  OF  G1(s)  WITH  T=2/3  &  1  SEC.  USING 

GAUSS 2AM 


FIGURE  7.4  NORMALIZED  NOISE  TRANSFERS  USING  PREDICTED 

PlCs  OF  G2(s)  WITH  8=1  AND  T=3/4  &  1  SEC. 

USING  GAUSS 2AM 
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sampling  rate  is  increased.  Therefore  the  noise  reduction 
for  the  system  containing  the  increased  sampling  rate  PIC 
in  place  of  the  normal  PIC  is  not  generally  as  large.  When 
the  value  of  the  pole  increases,  the  number  of  terms  in  the 
controller  has  less  effect  on  the  noise  transfer.  A  small 
change  in  the  sampling  rate  of  the  extended  controllers 
therefore  has  a  smaller  effect  than  at  the  lower  values  of 
the  pole. 

The  size  of  this  improvement  is  still  substantial. 
Therefore  for  applications  where  the  total  response  time  is 
more  important  than  the  noise  reduction  of  the  system,  this 
method  is  useful  to  consider. 

An  example  of  such  an  application  is  a  set-point 
system  where  the  level  of  the  input  may  change  fairly  often. 

In  this  case,  it  is  important  that  the  response  time  be  as 
low  as  possible.  Therefore  the  increased  sampling  rate  PlCs 
may  be  the  more  useful  type. 

This  method  may  be  extended  to  the  controllers  with 
two  or  more  extra  terms.  The  results  should  be  comparable  to 
those  obtained  for  the  PlCs. 

7-3  STRAIGHT  LINE  APPROXIMATIONS  TO  THE 
PREDICTED  PlCs 

In  some  situations,  it  is  useful  to  consider  methods 
to  approximate  the  predicted  PlCs.  By  using  an  approximation, 
a  more  rapid  method  of  calculation  is  available. 
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In  Chapter  3,  a  "straight  line"  approximation  is 
mentioned.  This  method  was  found  to  be  useful  for  particular 
plants.  The  approximation  minimizes  the  differences  between 
the  denominator  terms  of  the  controller  and  a  straight  line 
from  0  at  t  =  0  to  1.0  at  t  =  nT.  The  values  of  the  straight 
line  at  the  sampling  instants  are: 


n 


Therefore  the  "ideal"  errors  are: 


(7.1) 


(1  -  i),  <1  -  |),  •  •  •  ,  (1  -  Sjp)  (7.1a) 

The  differences  between  the  "ideal"  errors  and  the  error 
terms  of  the  controller  are: 

el>'  (1  -  I"  e2> . (1  "  ”  en) 

(7.2) 

Using  the  (7.2),  the  following  equation  may  be  written: 


Straight  Line  Difference  =  SLD  = 


(:L  "  n  -  el)2  +  (1  -  I  -  e2>2  + 


.  .  + 


/  n-1  \  2 

(1  -  —  "  en 
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For  each  plant,  the  error  terms  of  the  controller  e^r  e2' 

•  .  .  ,  e^  may  be  expressed  in  terms  of  e^.  Therefore  (7.3) 

may  be  rewritten  as: 

SLD  =  f(e±)  (7.4) 

The  value  of  e^  which  yields  the  absolute  minimum 
value  of  SLD  may  be  found  from  (7.4)  by  simple  differentiation 
as  in  Chapter  5.  Taking  the  derivative  of  (7.4)  with  respect 
to  e^  and  setting  the  result  to  zero  yields: 

aq-  tSLD1  =  a§7  [f(ei)J  =  0  <7-5> 

The  roots  of  (7.5)  correspond  to  the  local  maxima  and 

minima  of  SLD.  From  these  roots,  the  value  of  e^  which 

satisfies  the  deadbeat  requirements  and  corresponds  to  the 

minimum  of  SLD  is  chosen. 

The  other  coefficients  of  the  "straight  line" 

approximation  PIC  are  then  calculated  using  this  value  of  e^. 

In  FIG.  7.5  and  FIG.  7.6,  the  predicted  values  of  e^  are 

compared  to  the  "straight  line"  approximations  of  e^  for 

G-^(s)  and  (s)  respectively. 

From  FIG.  7.5  and  FIG.  7.6,  it  is  apparent  that  the 

approximation  is  most  accurate  for  <  1.  For  «  >  1,  the 

approximation  is  rather  inaccurate  and  therefore  not  useful. 

Although  the  results  are  not  shown,  the  approximation  is  also 

useful  for  a)  <  1.0  when  the  PIC  of  G „  (s)  is  used. 

n  4 
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a 

FIGURE  7.5  PREDICTED  e±  AND  STRAIGHT  LINE  APPROXIMATION 

e1  FOR  PlCs  OF  G1(s)  WITH  T=1  SEC. 

FIGURE  7.6  PREDICTED  e±  AND  STRAIGHT  LINE  APPROXIMATION 

e1  FOR  PlCs  OF  G2(s)  WITH  3=1  AND  T=1  SEC. 
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Similarly  this  approximation  may  be  extended  to  the 
P2Cs  where  the  approximate  ranges  of  usefulness  are  the 
same  as  for  the  PICs. 


CHAPTER  VIII 


GENERAL  CONCLUSIONS 

In  this  chapter,  some  general  conclusions  and 
observations  based  on  the  work  performed  for  this  thesis 
are  presented.  Then  some  extensions  of  this  work  and  the 
use  of  these  results  in  other  ways  are  discussed. 

8-1  CONCLUSIONS  AND  OBSERVATIONS 

A  design  criterion  for  the  calculation  of  non-minimum 
deadbeat  controllers  which  minimize  the  noise  transfer  of 
linear  systems  has  been  presented  in  this  thesis.  The  noise 
is  assumed  to  be  white  as  this  allows  the  application  of  the 
criterion  to  be  as  practical  as  possible. 

This  approach  to  the  noise  problem  is  novel  as  the 
additional  constraint  of  deadbeat  is  maintained  while  the  noise 
transfer  is  minimized  and  the  criterion  excludes  any  explicit 
reference  to  the  input  noise.  In  Section  3-5,  the  physical 
reasoning  behind  the  criterion  is  presented  as  part  of  the 
derivation  of  the  criterion. 

Two  types  of  extended  controllers  have  been 
considered: 

1.  Controllers  using  the  same  sampling  rate  as  the 
minimum  time  deadbeat  controllers. 

2.  Controllers  using  a  higher  sampling  rate  such 
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that  the  response  time  is  the  same  as  for  the 
minimum  time  deadbeat  controllers. 

For  the  first  case,  controllers  with  up  to  three 
extra  terms  were  considered.  As  stated  in  Section  5-7,  the 
noise  reduction  was  substantial  when  the  predicted  PIC  was 
used  in  place  of  the  minimum  time  controller.  From  Sections 
6-6  and  7-1,  the  extra  noise  reductions  obtained  by  using 
the  predicted  P2C  or  P3C  in  place  of  the  predicted  PIC  were 
not  as  large  as  the  initial  reduction. 

Therefore  when  the  possible  noise  reduction  is 
weighed  against  the  disadvantages  of  the  increased  response 
time  and  extra  controller  complexity,  the  optimum  practical 
extension  is  one  term.  However  for  applications  where  the 
absolute  noise  reduction  is  the  most  important  consideration, 
the  longer  extensions  will  be  useful. 

For  the  second  case,  only  the  PICs  for  G^(s)  and 
(s)  were  considered.  The  noise  reductions  obtained  using 
these  controllers  were  not  as  large  as  those  for  the  first 
type.  However  as  stated  in  Section  7-2,  the  reductions  were 
still  substantial.  Similar  results  should  be  obtainable  with 
other  plants.  A  similar  relationship  should  also  exist 
between  the  two  types  of  extended  controllers  when  other 
extensions  are  considered.  This  type  of  extended  controller 
is  most  useful  where  the  response  time  is  the  most  important 
consideration . 

The  PICs  predicted  by  the  design  criterion  are  a  good 
approximation  of  the  minimum  noise  transfer  controllers  for 


- 
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various  non-white  noise  sources.  This  was  verified  for  all 
plants  tested.  For  the  P2Cs  predicted  by  the  design  criterion, 
this  approximation  was  generally  not  as  good  as  for  the  PlCs. 

It  was  also  observed  that  some  predicted  PlCs  had  a 
lower  error  to  a  ramp  or  an  acceleration  input  than  did  the 
corresponding  minimum  time  deadbeat  controllers. 

Although  only  four  types  of  plants  were  tested,  the 
conclusions  reached  should  apply  to  any  linear  plant  in  a 
similar  system. 

8-2  EXTENSIONS  AND  DISCUSSION 

For  all  systems  tested,  the  noise  was  injected  before 
the  controller.  If  the  noise  is  injected  between  the  plant 
and  the  controller,  the  extended  controllers  with  an 
unchanged  sampling  rate  should  cause  only  a  slight  reduction 
in  the  noise  transfer  when  compared  to  the  noise  transfer 
obtained  using  the  minimum  time  controller.  Most  of  the 
noise  will  be  superimposed  directly  on  the  output  due  to  the 
slow  controller  response  time.  The  extended  controllers  with 
the  increased  sampling  rate  may  cause  slightly  more  improve¬ 
ment  because: 

1.  The  response  time  is  the  same  as  for  the 
minimum  time  deadbeat. 

2.  The  loop  gains  for  the  extended  controllers  are 
lower  than  for  the  minimum  time  controllers. 

3.  The  extended  controllers  can  sense  changes  at  the 
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output  more  quickly,  due  to  the  higher  sampling 
rate,  than  can  the  minimum  time  controllers. 

Although  the  calculation  of  a  design  criterion  would 
be  more  complicated,  the  method  should  be  applicable  to 
extended  deadbeat  controllers  designed  for  ramp  or  acceleration 
inputs  when  a  linear  plant  is  used.  It  should  still  be 
possible  to  minimize  the  "weighting  sequence"  of  an  extended 
deadbeat  controller  for  these  two  inputs. 

This  design  criterion  should  also  be  applicable  to 
systems  containing  a  linear  plant  and  a  linear  non-unity 
feedback  branch. 

If  an  otherwise  linear  system  contains  a  continuous 
saturable  type  of  non-linearity,  the  design  criterion  may  be 
of  some  use.  By  extending  the  controller,  the  gains  are 
lowered  and  thus  the  system  may  operate  in  the  linear  part 
of  the  non-linearity  during  a  greater  portion  of  the  time. 

As  long  as  the  non-linearity  operates  in  its  linear  region, 
the  noise  transfer  of  the  system  will  be  reduced.  For 
non-linearities  which  do  not  contain  linear  regions,  this 
design  criterion  would  not  be  useful. 
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APPENDIX  A 


HYBRID  COMPUTER  SIMULATION 

The  methods  used  to  simulate  the  plant  are  standard 
analog  techniques  and  therefore  will  not  be  illustrated  in 
this  appendix.  This  appendix  contains  only  the  details  of 
the  program  written  to  realize  the  digital  controller. 


A- 1  DIGITAL  CONTROLLER  REALIZATION 

The  digital  program  in  the  hybrid  simulation 
consisted  of  three  sections  which 

1.  Realized  the  controller 

2.  Realized  the  zero-order  hold 

3.  Calculated  the  noise  transfer  of  the  system  at 
the  end  of  the  test  run. 

The  realization  of  the  controller  function  was  hybrid 
The  coefficient  hQ  was  realized  on  the  analog  computer  while 
all  other  coefficients  appear  only  in  the  digital  program. 

A  digital  to  analog  converter  acted  as  the  zero-order  hold. 
The  diagram  of  the  realization  of  the  controller  and 
zero-order  hold  is  shown  in  FIG.  A.l. 

The  controller  was  realized  using  the  Direct  Digital 
Programming  discussed  in  Chapter  2.  It  was  shown  that  hQ 
operates  on  the  present  input  sample  to  the  controller.  All 
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other  coefficients  of  the  D(z)  operate  only  on  the  delayed 
information. 

Therefore  h^  was  realized  on  the  analog  computer  to 
make  the  length  of  the  sampling  window  as  short  as  possible 
for  the  equipment  used.  The  actual  length  of  the  sampling 
window  was  40  microseconds  in  the  hybrid  realization.  As 
shown  in  FIG.  A.l,  h^  BE1 (t)  was  summed  with  the  resultant 
products  of  the  delayed  information  and  the  other  controller 
coefficients . 

The  multiplications  of  the  controller  coefficients 
and  the  delayed  inputs  and  outputs  were  carried  out  between 
samples.  The  resulting  values  were  then  summed  and  the 
resultant  BE3(nT)  was  D/A  converted.  Then  at  the  sampling 
instants,  the  instantaneous  value  of  h^  BEl(t)  was  summed 
with  BE3(nT).  The  resultant  was  A/D  converted  and  then 
D/A  converted  immediately.  Both  the  preceding  and  following 
samplers  shown  with  the  controller  in  FIG.  1.1  were 
simulated  by  the  A/D  conversion.  The  zero-order  hold  was 
simulated  by  the  D/A  conversion.  The  length  of  the  sampling 
window  was  thus  set  by  this  A/D  and  D/A  sequence. 

To  relate  the  FOCAL  program  given  in  Section  A-2  to 
FIG.  A.l,  the  following  definitions  are  given: 


1. 

N 1  = 

n-1  where  n 

=  the  order  of  the  controller 

2. 

X  (i )  = 

-BEl* (t-iT) 

f  1  — •  1 ,  2  ^  •  •  •  ,  n 

3. 

Y(i)  = 

BE2  * (t-iT) 

j  1=1,  2 ,  •  •  • ,  n 

4. 

H(i)  = 

hi  ;  i  =  1, 

2  ,  •  •  •  ,  n 

. 
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5.  E(i)  =  e1  ;  i  =  1,  2 ,  .  .  n 
6  .  z  =  BE3 (nT) 

The  definitions  of  the  commands  used  in  FOCAL  are: 

1.  S  E  SET 

2.  T  E  TYPE 

3.  %  4.01  E  FORMAT  STATEMENT  (SAME  AS 
F  4.1  IN  FORTRAN) 

4.  FSQT  E  SQUARE  ROOT 

5.  FABS  E  ABSOLUTE  VALUE 

The  command  FNEW  is  used  in  FOCAL  to  perform  a 
variety  of  functions.  These  are  A/D  conversions,  D/A 
conversions  and  the  control  of  the  operating  state  of  the 
analog  computer  integrators. 

An  A/D  conversion  is  specified,  for  example,  by 

S  Y  (1)=FNEW (1,4,0) 

which  means  A/D  convert  the  signal  on  A/D  channel 
4  and  set  the  converted  value  equal  to  Y(l).  The  0  is  a 
dummy  variable  and  1  defines  an  A/D  conversion. 

A  D/A  conversion  is  specif ifed,  for  example,  by 

S  A=FNEW(2 ,3,Y (1) ) 


which  means  D/A  convert  the  value  Y (1)  onto  D/A  channel  3. 
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The  2  specifies  a  D/A  conversion  and  A  is  a  dummy  variable. 

The  command 

S  WAIT=FNEW(5 ,0,0,  ) 

is  used  to  synchronize  the  program  with  the  clock  of  the 
digital  computer.  This  clock  sets  the  length  of  the 
sampling  period  and  also  the  start  of  the  experimental  run. 
When  the  digital  computer  is  in  a  wait  loop,  the  only  way 
the  program  can  proceed  to  the  next  command  is  by  the 
reception  of  a  clock  pulse  which  occurs  at  the  sampling 
instants . 

The  FNEW  commands  which  control  the  operating  mode 
of  the  analog  integrators  are  self-explanatory. 

The  multiplications  of  the  controller  coefficients 
and  the  delayed  inputs  and  outputs  and  the  summing  of  the 
resultant  values  are  performed  by  statement  06.10. 

Statement  06.10  performs  both  the  multiplications 
and  the  summations  which  yields  BE3 (nT) .  The  shifting  and 
the  time  delaying  of  the  input  and  output  values  are  performed 
by  statements  06.50  to  06.80. 

When  the  analog  values  are  A/D  converted  in  statements 
07.15  and  07.20  to  allow  the  calculation  of  the  input  and 
output  root-mean-square  values,  the  divisor  2047  is  used. 

This  divisor  is  the  octal  value  in  the  A/D  converter  which 
corresponds  to  10  volts  on  the  analog  computer.  This 
division  is  necessary  to  yield  the  printed  values  of  NI  and 
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NO  in  volts.  These  are  the  root-mean-square  values  of  the 
input  and  output  noise  respectively. 


, 
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FIGURE  A. 1  HYBRID  REALIZATION  OF  THE  DIGITAL  CONTROLLER 


166 


A- 2  DIGITAL  PROGRAM  USED  IN  HYBRID  SIMULATION 


C-FOCAL. ,  1968 

01.10  S  H (1) = . . . ;  S  H (2 ) = . . . ;  ... 

01.20  S  E  ( 1 )  = .  .  .  ,  S  E  (2 )  =  . . . ;  ... 

01.70  S  Nl= . . . ;  S  P=100 ;  S  T=1 

02.10  S  M=0 ;  S  Z=0;  S  K=0;  S  El=FNEW (2,1,0) 
02.15  S  E3=FNEW (2,3,0) ;  S  WAIT=FNEW (5,0,0) 
02.20  S  OPERATE=FNEW (4,10,1) 

02.30  S  Y ( 1 ) =FNEW (1,4,0) ;  S  A=FNEW ( 2 , 3  ,  Y ( 1 ) ) 
02.32  S  X ( 1 ) =FNEW (1,1,0) 

02.35  FOR  R=1 , 1  ,K+1 ;  DO  6 . 1 
02.40  S  G=FNEW ( 2 , 1 ,  Z ) ;  S  Z=0 
02.45  IF  (K-l ) 6 . 8 , 6 . 7 
02.50  IF  (K-3 )  6 . 6 , 6 . 5 

03.10  S  M=M+ 1 ;  S  K=K+1;  IF  (P-M)7.1,7.1 
03.20  IF  (Nl-K) 4.1;  S  WAIT=FNEW (5,0,0) 

03.25  GO  TO  2.3 

04.10  S  WAIT=FNEW (5,0,0) 

04.20  IF  ( FABS ( FNEW (1,1,0) )-5)5. 1,5.1 
04.30  S  Y ( 1 ) =FNEW (1,4,0) ;  S  A=FNEW (2 , 3 ,Y (1) ) 
04.35  S  X ( 1 ) =FNEW (1,1,0) ;  S  K=0;  S  Z=0 
04.40  GO  TO  2.35 

05.10  S  E3=FNEW(2 ,3,0) ;  S  El^FNEW (2,1,0) 
05.20  GO  TO  3.1 
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06.10  S  B=H (R) *X (R) +E (R) *Y (R)+Z;  S  Z=B 
06.50  S  X  (5)  =X  (4  )  ;  S  Y(5)=Y(4) 

06.6  0  S  X  (4)  =X  (3)  ;  S  Y(4)=Y(3) 

06.70  S  X  ( 3 )  =X  (2  )  ;  S  Y(3)=Y(2) 

06.  80  S  X  (2)  =X  (1)  ;  S  Y(2)=Y(1);  GO  TO  3.1 

07.10  S  WAIT=FNEW(5/0,0) ;  S  HOLD=FNEW (4,10,4) 
07.15  S  NI=FSQT (10* (1/M*T) * (FNEW (1, 7, 0)/2047) ) 
07.20  S  NO=FSQT (10* (1/M*T)* (FNEW (1 , 8 , 0 ) /2047)  ) 
07.25  T  %12 . 06  ,  "  IN  "  ,NI,"  OUT  " ,NO,M 

07.30  T  %  12 . 06 ,  "  %  " ,  (N0*100)/Nl, !  I ! 

07.35  S  IC (INITIAL  CONDITION) =FNEW ( 4 , 10 , 2 ) 

08.10  QUIT 
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APPENDIX  B 


DIGITAL  COMPUTER  SIMULATION 
OF  THE  DIGITAL  CONTROLLERS 

The  digital  computer  program  used  to  simulate  the 
deadbeat  controllers  is  presented  in  this  appendix.  The 
program  sections  which  contain  the  variables  for  the  state 
equations  given  in  Chapter  1  and  the  calculation  algorithms 
for  the  minimum  time  and  extended  deadbeat  controllers  are 
not  listed.  These  sections  are  easily  programmed  directly 
from  the  lists  of  equations  given  in  the  text  of  the  thesis. 

B-l  DERIVATION  OF  PROGRAM 

This  program  incorporates  three  sections  which 

1.  Verify  the  deadbeat  response  of  the  system. 

2.  Subject  the  controllers  to  the  noise  sources 
defined  in  Chapter  4. 

3.  Calculate  the  noise  transfer  of  the  system 
containing  the  minimum  time  deadbeat  controller 
and  the  normalized  noise  transfers  of  the 
extended  deadbeat  controllers. 

These  sections  are  numbered  in  the  program  listing. 

The  systems  are  realized  by  combining  the  discrete 
state  equations  with  the  direct  digital  realization  of  the 
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digital  controller.  In  this  program,  AH (N )  and  AE (N)  are 
used  to  represent  and  e^  respectively. 

B-2  DIGITAL  PROGRAM 


100  CONTINUE 

C  AIN 2  AND  AOUT2  ARE  USED  TO  CALCULATE 

C  INPUT  AND  OUTPUT  MEAN-SQUARE  NOISE 

AIN  2=0 . 0 
AOUT  2=0.0 

C  SET  XI (1)  TO  XN (1)  TO  ZERO 

X1(1)=0.0 

• 

XN ( 1 ) =  0 . 0 

C  INITIALIZE  OTHER  COEFFICIENTS 

ICYCLE=N 
K=1 
JJ=1 

COUNT=l. 0 
XOUT1  (1) =0 . 0 

C  CONTROLLER  REALIZATION  BY  DIRECT  DIGITAL 

C  PROGRAMMING 

C  CALCULATE  INPUT  TO  CONTROLLER , BE 1 (JJ) 

101  BEl ( JJ) =1 . 0-XOUT1 ( JJ) 

BEA1=0 . 0 

BEB1=0 . 0 

DO  560  JK=1 ,  K 

BEA ( JK) =AH (JK) *BE1 (JJ-JK+1) 

BEAl=BEAl+BEA ( JK) 

560  CONTINUE 

IF (K.LT. 2)  GO  TO  562 
DO  561  JL=2 ,  K 

BEB (JL) =AE ( JL) *BE2 ( JJ-JL+1) 

BEBl=BEBl+BEB ( JL) 
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561  CONTINUE 

C  CALCULATE  OUTPUT  OF  CONTROLLER,  BE 2 (JJ) 

562  BE2 ( JJ) =BEA1+BEB1 

C  STATE  EQNS .  TO  CHECK  DEADBEAT  (SECTION  1) 

Xl (K+l) = .  .  . 

• 

XN (K+l) = .  .  . 

JJ=JJ+1 

XOUT1 (JJ)=X1 (K+l) 

IF  (COUNT  . GE .7.0)  GO  TO105 
COUNT=COUNT+l . 0 
IF (K.GE.ICYCLE)  GO  TO102 
K=K+1 
GO  TO101 

102  XI  (1)=K1  (K+l) 

• 

XN ( 1 ) =XN  (K+l) 

K=1 

GO  TO101 
105  X1(1)=0.0 

• 

XN  (1) =0 . 0 

K=1 

J=1 

XOUT ( 1 ) =0 . 0 

C  CALCULATE  INPUT  TO  CONTROLLER,  BEl(J) 

C  WHERE  YFLl(J)  IS  THE  NOISE  SOURCE  VALUE 

111  BE1  (J)=YFL1 (J) -XOUT (J) 

BEB1=0 . 0 

BEA1=0 . 0 

DO  550  JM=1 , K 

BEA  ( JM)  =AH  ( JM)  *BE1  ( J- JM+1) 

BEA1=BEA1+BEA  ( JM) 

550  CONTINUE 

IF  (K.LT.2)  GO  T0552 
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DO  551  JN=2,K 

BEB (JN)=AE (JN) *BE2 (J-JN+1) 

BEB1=BEB1+BEB (JN) 

551  CONTINUE 

C  CALCULATE  OUTPUT  OF  CONTROLLER,  BE 2 (J) 

552  BE2  (J)=BEA1-BEB1 

C  STATE  EQNS .  TO  TEST  SYSTEM  IN 

C  PRESENCE  OF  NOISE  INPUT  (SECTION  2) 

Xl  (K+l )  =  .  .  . 

• 

XN (K+l ) = .  .  . 

J=J+1 

XOUT (J)=X1 (K+l) 

C  NRUN=256  OR  1024,  ICYCLE=N 

IF (J.GE. (NRUN+1) )  GO  T0113 
IF (K . GE . I CYCLE)  GO  T0112 
K=K+1 
GO  TOlll 

112  XI (1)=X1 (K+l) 

• 

XN ( 1 ) =XN (K+l) 

K=  1 

GO  TOlll 

C  CALCULATE  THE  SUM  OF  SQUARES  OF  THE 

C  NOISE  IN  AND  OUT  (SECTION  3) 

113  DO  114  1=1,  NRUN 
AIN2=AIN2+YFL1  (I) **2 
AOUT2=AOUT2+ (XOUT  (1  +  1) ) **2 

114  CONTINUE 

C  CALCULATE  MEAN-SQUARE  TRANSFER  RATIO 

RMSNSE=AOUT 2/AIN  2 

C  IF  CHECK=0 . 0 ,  RETURN  TO  MINIMUM  TIME 

C  CONTROLLER  PROGRAM 

IF (CHECK. EQ. 0 . 0)  GO  TO150 
C  CALCULATE  NORMALIZED  NOISE  TRANSFER 
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C  FOR  EXTENDED  CONTROLLER  WHERE  RMSDBT 

C  IS  THE  NOISE  TRANSFER  OF  MINIMUM 

C  CONTROLLER 

RMS=RMSNSE/RMSDBT 

C  RETURN  TO  EXTENDED  CONTROLLER  PROGRAM 

GO  TO170 
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